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I.  Introduction 

This  report  contains  user  instructions,  a  listing  and 
documentation  for  a  microcomputer  BASIC  program  that  can  be  used 
to  compute  an  estimate  of  the  probability  that  a  magnetometer 
based  detection  system  such  as  the  AN/ASQ-81  will  detect  a 
submarine  magnetic  dipole  field  during  an  encounter. 

The  program  generates  detection  probabilities  based  on  two 
encounter  models.  In  the  first  encounter  model,  the  detection 
system  uses  a  crosscorrelation  detector.  In  the  second  encounter 
model,  the  detection  system  uses  a  square  law  detector.  Relative 
to  operationally  realizable  values,  probabilities  based  on  the 
first  model  represent  upper  bounds  and  those  based  on  the  second 
model  represent  lower  bounds.  For  both  encounter  models,  the 
signal  is  proportional  to  the  square  of  the  magnitude  of  the 
projection  of  the  dipole  field  on  the  earth  magnetic  field  and 
the  magnetic  noise  does  not  change  with  changes  in  the  position 
of  the  magnetometer.  Also,  for  both  encounter  models,  an 
encounter  is  a  straight  line  encounter  with  constant  vertical 
separation. 

The  encounter  models  can  be  interpreted  as  models  of  a 
magnetic  anomaly  detection  (MAD)  system  on  an  aircraft  that  is 
moving  with  constant  course,  speed  and  altitude  in  an  encounter 
with  a  submarine  moving  with  constant  course,  speed  and  depth. 

Or,  they  can  be  interpreted  as  models  of  a  stationary  magnetic 
anomaly  detection  system  in  an  encounter  with  a  submarine  moving 
with  constant  course,  speed  and  depth. 
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The  program  parameters  include:  encounter  latitude  and 
longitude,  submarine  induced  magnetic  moments,  submarine 
permanent  magnetic  moments,  submarine  course,  speed  and  altitude, 
magnetometer  course,  speed  and  depth,  encounter  lateral  range 
(the  horizontal  range  at  the  closest  point  of  approach  in  a 
straight  line  encounter)  and  false  alarm  rate.  (In  the  program, 
a  false  alarm  is  the  event  that  the  detection  system  classifies 
noise  as  a  dipole  signal.) 


II.  Program  Usar  Instructions 

As  listed  in  Appendix  10 ,  the  program  can  be  run  under  a 
BASIC  language  that  is  compatible  with  IBM  PC  BASIC.  If  the 
listing  is  used  to  enter  the  program  through  the  keyboard,  then 
the  program  should  be  saved  with  the  name  MAD. BAS  or  the  value  of 
N$  on  line  40  should  be  changed  to  the  file  name  under  which  it 
is  saved. 

The  program  contains  user  instructions  in  the  form  of  query 
and  parameter  limitation  messages.  As  an  example  of  the  former, 
after  starting  the  program  under  BASIC,  the  following  message 
should  appear: 

Magnetic  Anomaly  Detection  (MAD)  Lateral  Range  Function 

generate  or  print  a  program  data  file  (g/p)? 

By  entering  p,  data  can  be  printed  from  a  program  data  file 
that  was  generated  by  the  program.  By  entering  g,  a  program 
data  file  can  be  generated  for  a  set  of  user  specified 
conditions.  With  either  response,  a  sequence  of  additional 
queries  is  displayed.  These  queries  require  either  an  indicated 
response  or  a  parameter  value  as  the  input.  If  the  initial 
response  is  g,  the  sequence  includes  queries  whose  responses 
determine  whether  or  not  an  auxiliary  data  file  will  be  generated 
that  can  be  used  for  future  input  of  magnetic,  processing  or 
kinematic  data.  In  particular,  the  first  query  in  the  sequence 
gives  the  option  of  using  a  combined  magnetic,  processing  and 
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kinematic  data  file.  If  the  response  indicates  that  it  should  be 
used  and  the  file  is  available,  the  parameter  values  that  remain 
to  be  entered  in  order  to  generate  a  program  data  file  are  the 
following:  the  false  alarm  rate,  the  magnetic  noise,  the  maximum 
encounter  lateral  range  and  the  lateral  range  step.  The  combined 
file  should  be  used  only  if  the  effect  of  varying  just  one  or 
more  of  these  parameters  is  desired.  If  the  response  indicates 
that  the  file  will  not  be  used,  queries  concerning  magnetic, 
kinematic  and  processing  parameter  values  are  displayed. 

After  all  of  the  program  parameter  values  have  been  entered, 
a  query  is  displayed  giving  the  option  of  generating  the  combined 
file.  Then  a  query  giving  the  option  of  generating  a  program 
data  file,  a  query  giving  the  option  of  printing  the  encounter 
parameter  values  and  a  query  giving  the  option  of  printing 
lateral  range  function  values  are  displayed.  The  lateral  range 
function  values  are  the  encounter  detection  probabilities  indexed 
by  lateral  range.  The  parameters  maximum  lateral  range  and 
lateral  range  step  determine  the  index  lateral  ranges  of  the 
encounters  for  which  probabilities  of  detection  are  computed. 

The  program  generates  magnetic  signal  values  that  correspond 
to  points  in  time  during  an  encounter.  Following  the  lateral 
range  function  query,  a  query  is  displayed  that  gives  the  option 
of  printing  the  magnetic  signal  values  for  an  encounter.  If  the 
option  is  exercised,  the  option  is  repeated.  When  the  option  is 
not  exercised,  a  query  is  displayed  that  gives  the  option  of 
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generating  or  printing  a  new  program  data  file.  If  this  option 
is  not  exercised,  the  program  ends. 

Some  suggested  guides  for  determining  parameter  values  can 
be  found  in  Section  IV  of  this  report. 
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III.  Encounter  Model  Limitations 

In  the  two  encounter  models,  the  signal  is  the  unfiltered 
signal  that  is  generated  by  a  magnetic  dipole  that  moves  relative 
to  a  magnetometer.  Depending  on  the  input  filtering,  describing 
a  ship  magnetic  anomaly  as  a  dipole  anomaly  should  not  be  a 
significant  limitation  for  an  encounter  slant  range  at  the 
closest  point  of  approach  (CPA)  that  is  greater  than  one  hull 
length.  The  detection  decision  is  based  on  samples  from  a  single 
time  interval  (window)  that  is  centered  on  the  CPA.  The  length 
of  the  sampling  interval  and  the  sampling  rate  are  parameter 
values  that  are  inputs  to  the  program.  In  terms  of  signal-to- 
noise  ratio,  there  is  an  optimum  sampling  interval  length 
(integration  time)  and  sample  rate.  Although  dipole  signal 
energy  is  not  symmetrically  distributed  about  the  CPA  time,  for  a 
given  sampling  interval,  the  difference  between  the  signal  energy 
for  the  optimum  interval  location  and  the  CPA  centered  location 
should  not  be  significant  in  most  cases. 

The  encounter  model  magnetometer  noise  samples  are 
determined  by  a  gaussian  random  process.  They  are  values  of 
independent  identically  distributed  normal  random  variables.  The 
standard  deviation  of  these  random  variables  is  referred  to  as 
the  magnetic  noise  and  the  variance  is  the  noise  in  the  sense  of 
the  signal-to-noise  ratio.  The  degree  of  correspondence  between 
this  process  and  operational  noise  depends  on  the  nature  of  the 
dominant  operational  noise  sources  and  on  the  magnetometer  input 
filtering  (noise  whitening). 
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In  the  encounter  models,  the  intervals  are  adjacent  but  not 
overlapping  and  a  detection  statistic  corresponds  to  each  sample 
interval  and  its  value  is  determined  by  the  sample  values.  If 
the  value  of  the  detection  statistic  for  an  interval  equals  or 
exceeds  a  threshold  value,  a  detection  is  indicated.  The 
threshold  value  is  determined  by  the  false  alarm  probability 
which  in  turn  is  determined  by  the  false  alarm  rate  and  the 
sampling  interval  length. 

A  detector  that  used  a  moving  sample  interval  that  was 
generated  by  replacing  the  oldest  sample  by  the  newest  one  would 
correspond  more  closely  to  the  detector  in  an  operational 
detection  system.  In  a  model  of  the  detector,  since  the  sample 
windows  overlap,  the  detection  statistics  would  represent  a 
sequence  of  random  variables  that  were  correlated  over  an 
interval  equal  to  the  width  of  the  sample  interval.  Because  of 
this  dependence,  it  seems  unlikely  that  the  results  that  would  be 
obtained  with  an  encounter  model  based  on  an  overlapping  interval 
detector  would  differ  significantly  from  those  obtained  with  an 
encounter  model  based  on  a  nonoverlapping  interval  detector. 

For  encounter  lengths  of  the  order  of  a  few  nautical  miles 
or  less,  the  straight  line  encounter  condition  should  not  be  a 
significant  limitation.  In  particular,  this  should  be  the  case 
for  a  fixed  magnetometer  since,  for  a  submarine  (or  surface  ship) 
target,  vertical  separation  and  course  changes  should  be  le»s 
likely  to  occur. 
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Other  models  are  available  that  can  be  used  as  the  basis  for 
computing  an  estimate  of  the  probability  that  a  magnetic  anomaly 
detection  system  will  detect  a  submarine  during  an  encounter. 

For  example,  one  that  is  described  in  Appendix  8  can  be  used  to 
determine  the  slant  ranges  of  straight  line  encounters  for  which 
the  detection  probability  is  equal  to  a  specified  value.  The 
parameter  values  that  are  required  to  do  this  are  an  average 
submarine  dipole  moment,  a  detection  system  capability  factor  and 
a  noise  factor.  Values  for  these  parameters  can  be  determined 
from  operational  data.  However,  the  values  are  specific  to 
averages  over  a  particular  set  of  encounter  conditions.  An 
advantage  of  the  two  encounter  models  relative  to  this  model  is 
their  adaptability  to  different  magnetic,  processing  and 
kinematic  conditions. 
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IV.  Parameter  Values 

The  magnetic  parameter  value  queries  are  generally 
explanatory  with  regard  to  the  value  that  should  be  entered. 

This  is  also  true  of  the  kinematic  parameter  value  queries. 
However,  there  is  some  ambiguity  with  respect  the  processing 
parameter  value  queries  and  the  noise  parameter  value  query.  To 
reduce  this  ambiguity,  a  brief  discussion  of  the  common 
characteristics  of  the  two  encounter  model  processing  parameters 
and  the  noise  parameter  is  given  below.  This  is  followed  by  some 
guidelines  for  choosing  these  parameter  values. 

In  both  encounter  models,  a  decision  is  made  at  the  end  of 
each  sampling  interval.  The  decision  is  either  noise  energy  was 
present  during  the  interval  or  noise  energy  plus  signal  energy 
was  present  during  the  interval.  The  sample  intervals  are 
adjacent,  equal  width,  nonoverlapping  time  intervals.  The  number 
of  samples  that  are  input  in  a  sample  interval  is  determined  by 
the  sampling  rate  and  the  interval  length. 

The  program  default  choice  for  the  sampling  rate  is  2 -MAXF 
where  MAXF  is  a  parameter  that  is  labeled  the  maximum  magnetic 
signal  frequency.  This  sample  rate  is  the  Nyquist  rate  for  an 
ideal  low  pass  filter.  However,  tha  signal  in  a  sample  interval 
that  is  computed  by  the  program  represents  an  unfiltered  dipole 
signal.  This  is  a  reasonable  approximation  if  the  signal  energy 
that  is  associated  with  signal  components  greater  than  MAXF  is 
relatively  small.  As  discussed  below,  the  noise  energy  in  a 
sample  interval  should  be  considered  to  be  proportional  to  MAXF 
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in  order  to  be  consistent  with  the  encounter  models.  Ideally,  a 
default  choice  for  MAXF  would  make  the  ratio  of  the  signal  energy 
to  the  noise  energy  a  maximum  for  the  sampling  interval  of  an 
encounter.  The  program  default  choice  for  MAXF  is  2  MAXVM/MINRO 
where  MAXVM  is  a  user  estimated  maximum  encounter  relative  speed 
converted  from  knots  to  meters  per  second  and  MINRO  is  a  user 
estimated  maximum  slant  range  at  CPA  in  meters  in  terms  of  a  just 
detectable  target.  (A  precise  definition  of  a  just  detectable 
target  can  be  made  in  terms  of  a  specified  detection  probability, 
false  alarm  probability  and  target  dipole  moment.)  The  default 
choice  for  MA  IF  is  consistent  with  the  observation  in  Reference  1 
that  if  an  optimum  value  for  MAXF  is  determined  for  a  minimum 
dipole  moment  target,  then  no  significant  increase  in  MAXF  is 
required  in  order  to  maintain  a  required  detection  probability  if 
the  encounter  lateral  range  is  decreased  even  though  the  signal 
energy  spectrum  is  shifted  to  higher  frequencies. 

Because  of  the  detection  statistics  and  the  gaussian  noise 
model,  if  there  were  no  penalty  for  decision  delay,  a  sample 
interval  length  for  a  signal  should  be  chosen  equal  to  the  signal 
duration,  since  this  would  make  the  detection  probability  for  an 
encounter  a  maximum.  In  the  program,  a  default  choice  for  the 
sample  interval  length  (integration  time)  is  2 -MAXRO/MINVM  where 
MINVM  is  a  user  estimated  minimum  encounter  relative  speed 
converted  from  knots  to  meters  per  second  and  MAXRO  is  a  user 
estimated  maximum  slant  range  at  CPA  in  meters  for  a  detectable 
encounter  (in  terms  of  a  specific  detection  probability  and  false 
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alarm  probability)  with  a  user  specified  maximum  dipole  moment 
target.  This  choice  might  be  considered  a  balance  between 
minimizing  decision  delay  and  maximizing  detection  probability. 

In  the  encounter  models,  the  sample  intervals  are  located  so  that 
the  CPA  time  is  at  the  center  of  a  sample  interval  and  this 
sample  interval  is  the  only  one  that  contains  signal  energy. 

This  characteristic  is  consistent  with  the  default  choice  for  the 
sample  interval.  For  a  given  sample  interval  length,  although  in 
general  the  interval  location  is  not  the  optimum  one  in  terms  of 
signal  energy,  it  should  be  approximately  so  in  most  cases. 

The  program  noise  parameter  is  SIG.  It  represents  the 
standard  deviation  a  associated  with  the  magnetic  noise  process 
of  the  two  encounter  models.  In  terms  of  the  ideal  low  pass 
filter  implied  by  the  encounter  models,  its  square  should  be 
equal  to  MAXF-(SIGO)2  where  SIGO  is  the  magnetic  noise  process 
constant  power  spectral  density.  The  program  does  not  enforce 
this  relation.  Therefore,  in  using  the  program,  the  implied 
relation  between  the  two  input  parameters:  magnetic  noise  and 
maximum  magnetic  signal  frequency  should  be  kept  in  mind.  If  ah 
average  value  of  the  peak-to-peak  magnetic  noise  for  an  encounter 
can  be  estimated,  for  example  from  a  magnetometer  trace,  then  the 
value  for  SIG  should  be  chosen  so  that  the  estimate  is  4  to  6 
times  this  value. 
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Appendix  1.  The  Detection  Statistics 

In  this  appendix,  y,,  y2,  •  •  ■  ,  yM  are  sequential  values 
(voltages)  representing  the  sample  values  in  a  sample  interval. 
They  are  the  input  to  a  magnetometer's  detector.  With  these  M 
sample  values,  the  detector  computes  the  value  of  a  detection 
statistic.  This  value  is  represented  by  x  and  the  detection 
decision  corresponding  to  the  sample  interval  is  determined  by 
the  decision  rule:  If  x  >  x‘,  then  the  input  during  the  sample 
interval  was  noise  plus  signal,  otherwise,  the  input  was  noise. 

For  both  encounter  models,  the  detection  probability  and  the 
false  alarm  probability  are  decreasing  functions  of  x*  and  the 
relation  is  one-to-one  in  both  cases.  In  the  program,  the  false 
alarm  probability  pF  is  used  to  determine  a  unique  value  of 
the  threshold  x*.  This  value  is  then  used  to  determine  a  unique 
value  of  the  detection  probability  pD. 

In  the  program,  pF  is  found  using  the  relation  pF  =  R-6t 
where  R  is  the  false  alarm  rate  in  false  alarms  per  second 
and  6t  is  the  sample  interval  length  in  seconds.  This  relation 
is  based  on  the  following  argument:  With  no  signal  energy  in  a 
sample  interval,  y,  =  n,,  ys  =  n2,  •  •  •  ,  y  =  nM  where  n,,  n2,  ••• 

,  nM  are  noise  values  (voltages)  input  to  a  magnetometer's 
detector  that  are  values  of  independent  normal  (gaussian)  random 
variables,  each  with  mean  zero  and  standard  deviation  a. 

Because  of  this,  in  the  encounter  models,  values  of  x  for 
different  sample  intervals  are  the  values  of  independent  random 
variables  that  determine  two  outcomes:  x  >  x’  or  x  <  x‘. 
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Therefore,  pF  is  the  same  for  each  sample  interval  and,  in 
terms  of  these  outcomes,  the  sequence  is  a  series  of  repeated 
independent  Bernoulli  trials  and  the  expected  number  of  trials 
between  false  alarms  is  l/pF.  Since  the  time  between  trials  is 
£t,  the  expected  number  of  seconds  between  false  alarms  is  6t/pF 
and  the  expected  number  of  false  alarms  per  second  is  equal  to 
pF/$t.  So,  R  =  pf/St. 

The  determination  of  x*  depends  on  the  encounter  model 
statistic.  For  both  encounter  models,  when  there  is  a  signal, 
y,  -  n,  +  s,,  y,  ■  n,  +  y„  •  •  •  ,  yM  -  nM  +  sM  where  st,  s2,  •  •  •  , 
sM  are  signal  values  (voltages)  input  to  a  magnetometer's 
detector.  The  models  imply  that  the  signal  values  s  *  c-Hs 
where  c  is  a  constant  whose  value  is  determined  by  the 
characteristics  of  the  encounter  magnetometer  and  where  the  Hs 
are  dipole  magnetic  signal  intensities.  The  models  also  imply 
that  the  noise  values  n  are  determined  by  a  gaussian  stochastic 
processes  characterized  by  a  standard  deviation  a  and  that 
n  =  c-Hn  +  n'  where  the  HN  are  magnetic  noise  intensities  and 
the  n'  are  magnetometer  instrument  noise  values  (voltages) .  In 
the  program,  the  magnitude  of  c  is  1.  Since,  for  both  models, 
pD  depends  only  on  the  ratio  of  signal  energy  to  noise  energy 
for  a  sample  interval,  this  is  a  satisfactory  choice  for  the 
program  if  instrument  noise  process  is  assumed  to  be  independent 
of  the  magnetic  noise  process  and  to  be  expressed  in  terms  of  an 
equivalent  magnetic  noise  H  =  (1/c)  n'. 
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For  both  encounter  models,  the  signal  (the  average  signal 
power)  S  =  (l/M)  *E  sK2  and  the  noise  (the  expected  value  of  the 
average  noise  power)  N  =  a*  so  that  the  signal-to-noise  ratio 
is  (1/M)  -E  s K*/a2  where  the  sum  index  K  ■  1,  2,  •••  ,  M. 

The  Crosscorrelation  Detector  Statistics  The  statistic  for  the 
first  encounter  model  is  a  crosscorrelation  detector  statistic 
that  is  defined  by  the  sum 

x  *  E  yK-sK 

where  the  summation  index  K  =  1,  2,  •••  ,  M  and  the  sum  is  over 

the  values  corresponding  to  a  sample  interval.  For  the  first 
encounter  model,  the  characteristics  of  both  the  noise  and  the 
signal  are  required  in  order  to  determine  encounter  detection 
probabilities.  In  particular,  the  signal  values  for  an  encounter 
are  in  the  memory  of  the  detector  prior  to  the  encounter.  For 
the  encounter  conditions  and  a  specified  false  alarm  probability, 
the  statistic  is  optimum  in  the  sense  that  the  encounter 
detection  probability  for  this  statistic  is  at  least  equal  to 
that  for  any  other  statistic.  Because  of  these  considerations, 
encounter  probabilities  based  on  the  crosscorrelation  statistic 
can  be  considered  to  represent  upper  bounds  on  detection 
performance  against  dipole  targets  for  magnetometers  of  the  type 
described  by  the  models. 

For  a  sample  interval  without  signal  energy,  x  is  the  value 
of  a  normal  random  variable  with  a  mean  jix  =  0  and  a  variance 
ax2  =  o*  •  E  sK2  where  a  is  the  standard  deviation  associated  with 
the  noise  process  and  the  sum  index  K  =  1,  2,  • • •  ,  M  and  the 


14 


sum  is  over  the  values  corresponding  to  the  sample  interval. 

This  implies  that 

Pf  -  1  -  P(xVCTx) 

where  P(z)  is  the  standard  normal  cumulative  distribution 
function.  This  relation  is  the  basis  for  determining  the 
threshold  value  x*. 

For  the  sample  interval  with  signal  energy,  x  is  the  value 
of  a  normal  random  variable  with  a  mean  =  Z  sK*  where  the  sum 
index  K  -  1,  2,  • • •  ,  M  and  the  sum  is  over  the  values 

corresponding  to  the  sample  interval.  This  implies  that 

p0  *  1  -  P(v*  -  dVi) 

where  v*  =  x‘/crx  and  d  =  Z  sK2/ff*  “  M*S/N.  This  relation  is  the 
basis  for  determining  encounter  detection  probabilities  for  the 
first  encounter  model.  The  relation  implies  that  for  a  specified 
false  alarm  probability  pf  the  detection  probability  pD  is  an 
increasing  function  of  the  signal  to-noise  ratio  S/N. 

The  Square  Law  Detector  Statistic:  The  statistic  for  the  second 
encounter  model  is  a  square  law  (energy)  detector  statistic  that 
is  defined  by  the  sum 

x  -  Z  yK* 

where  the  sum  index  K  =  1 ,  2 ,  • • •  ,  M  and  the  sum  is  over  the 

values  corresponding  to  the  sample  interval.  For  the  second 
encounter  model,  only  the  characteristics  of  the  noise  are 
required  to  determine  encounter  detection  probabilities. 
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For  a  sample  interval  without  signal  energy,  x/o*  is  the 
value  of  a  chi-square  random  variable  with  M  degrees  of 
freedom.  This  implies  that 

Pf  ■  1  -  P(x‘/a*|M) 

where  P(x‘/az|M)  is  the  chi-square  cumulative  distribution 
function  for  a  chi-square  random  variable  with  H  degrees  of 
freedom  and  where  a  is  the  standard  deviation  associated  with 
the  noise  process.  This  relation  is  the  basis  for  determining 
the  threshold  value  x‘. 

For  the  sample  interval  with  signal  energy,  x/a*  is  the 
value  of  a  noncentral  chi-square  random  variable  with  M  degrees 
of  freedom  and  noncentral  parameter  £  sK*/a2  where  the  sum  index 
K  =  1 ,  2 ,  • ■ •  ,  M  and  the  sum  is  over  the  values  corresponding 

to  the  sample  interval.  This  implies  that 

Po  =  1  ~  P(x'/o2|M,  e  sKVa2) 

where  P(x*/ct7|m,  E  sKVa?)  is  the  cumulative  distribution 
function  for  a  noncentral  chi-square  random  variable  with  M 
degrees  of  freedom  and  noncentral  parameter  E  s^/a1-  M-S/N. 

This  relation  is  the  basis  for  determining  encounter  detection 
probabilities  for  the  second  encounter  model.  The  relation 
implies  that  for  a  specified  false  alarm  probability  pF,  the 
detection  probability  pD  is  an  increasing  function  of  the 
signal-to-noise  ratio  S/N.  This  is  made  more  evident  by  the 
following  relation: 

P(x‘/a’|M,  E  s’/o1)  =  E{  (aJ/J!)  -exp(-a)  -P[xVaz|  (M  +  2  •  J)  ]  } 
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where  the  parameter  a  =  (1/2)  *2  s £/o%  ■»  (M/2)  *(S/N),  the  sum 


index  J  =  0 ,  1 ,  2 ,  • • •  and  the  sum  index 

Note,  since  P(x*/a*|M)  <  P[x‘/o2|  (M  +  2’J)] 

,  the  relation  p0  £  pF  is  satisfied. 


K  *  1,  2,  •  •  •  ,  M. 

for  J  =  0,  1,  2, 
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Appendix  2.  Program  Probability  Calculations 

The  program  evaluates  the  cumulative  and  inverse  cumulative 
distribution  functions  using  approximations  described  in 
Reference  2.  These  approximations  are  listed  below. 

Ths  Standard  Normal  cumulative  Distribution  Function 
Approximation : 

P(z)  =  1  -  s-t*(b,  +  t-(b2  +  t*(bs  +  t  •  (b4  +  t -b5) ) ) ) 
where  s  =  (1/(2  -it)'1*)  -exp(-z*/2)  and  t  *  1/(1  +  b0*z).  And 
where 

b0  -  .2316419,  b,  =  .319381530,  b2  =  -.356563782, 
b4  =  -1.821255978,  b,  =  1.330274429, 
and  z  >  0.  For  z  <  0,  P(z)  =  1  -  P(|z|). 

The  Inverse  standard  Normal  Cumulative  Distribution  Function 
Approximation: 

z(P)  =  t  -  ( c0  +  t-(c,  +  t  ■  c2) )  /  ( 1  +  t-(d,  +  t-(d2  +  t-d3))) 
where  t  =  (ln(l/Q2) ) 1/2  and  Q  =  1  -  P.  And  where 
c0  =  2.515517,  c,  =  .802853,  C2  =  .010328, 
d,  =  1.432788,  d2  -  .189269,  d3  =  .001308, 
and  .5  <  P  <  1.  For  0  <  P  <  .5,  z(P)  =  -z(l  -  P) . 

The  inverse  Chi-Square  cumulative  Distribution  Function 
Approximation: 

v(P|M)  =  M  •  [  1  -  2/(9  -M)  +  Z  •  (2/(9 -M)  )1/2]3 
where  P(z)  =  P(v|M).  In  the  program,  the  inverse  standard 
normal  cumulative  distribution  function  approximation  is  used  to 
determine  z. 
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The  Moncentral  Chi-Square  cumulative  Distribution  Function 
Approximation : 

P(w|M,E  sk8/c!)  =  P(z) 

where  z  *  [2-w/(l  +  b)]1/a  -  [2>a/(l  +  b)  -  l]v*  with 

a  «  M  +  2  sj/o*,  b  =  (2  sKVc2)/(M  +  2  sKa/a*)  and  the  sum  index 
K  =  l,  2 ,  •••  ,  M.  In  the  program,  the  standard  normal 
cumulative  distribution  function  approximation  is  used  to 
determine  P(z) . 
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Appendix  3.  The  Magnetic  signal 

The  encounter  models  are  defined  by  the  following 
conditions:  A  submarine  magnetic  anomaly  field  is  a  magnetic 
dipole  field  that  is  superimposed  on  an  earth  magnetic  field  that 
is  constant  over  an  encounter  region.  A  magnetometer  magnetic 
signal  value  is  the  magnitude  of  the  projection  of  a  dipole 
magnetic  field  on  the  earth  magnetic  field  at  the  location  of  the 
magnetometer.  The  basis  for  determining  magnetic  signal  values 
in  the  program  is  an  expression  that  can  be  developed  as 
follows:  In  the  right  handed  rectangular  coordinate  system  that 

is  shown  in  Figure  1,  a  magnetic  dipole  is  at  the  origin,  the  xy- 
plane  is  the  horizontal  plane  at  a  representative  point  in  an 
encounter  region,  the  positive  y-axis  is  in  the  direction  of 
magnetic  north,  the  positive  x-axis  is  in  the  direction  of 
magnetic  east  and  the  positive  z-axis  is  positive  upward.  In 
this  rectangular  coordinate  system,  the  constant  earth  magnetic 
field  can  be  expressed  by  HE  =  HE-(j-cos  *  +  It  sin  #)  where  * 
is  a  magnetic  dip  angle  and  HE  is  a  magnetic  field  magnitude 
that  characterizes  the  earth  field  in  an  encounter  region. 

In  a  spherical  coordinate  system  with  the  origin  at  the 
magnetic  dipole  and  the  polar  axis  in  the  direction  of  the  dipole 
moment,  Hd  =  (c -p/r3)  •  (2 -r -cos  6  +  0-sin  9)  is  the  magnetic 
field  of  the  dipole  at  a  point  whose  spherical  coordinates  are 
(r,r,9).  In  this  expression  p  is  the  magnitude  of  the  dipole 
moment  and  c  is  a  constant  whose  value  is  determined  by  the 
choice  of  units.  The  magnetic  signal  for  a  magnetometer  that  is 


20 


Figure  1.  A  unit  vector  in  the  direction  of  the  earth  field  and 
unit  vector  in  the  direction  of  the  dipole  field  of  a 
dipole  at  the  origin  are  shown  at  a  point  that  is  a 
distance  r  from  the  dipole.  The  unit  vector  h  is 
in  the  direction  of  the  horizontal  component  of  the 
dipole  moment. 


deaoribed  by  the  encounter  models  i«  H,  ■  when  the 

magnetometer  in  At  the  point  (r,r,b)  and  H,  een  be  expressed  in 
terms  of  the  rectangular  ooordinete  ayetsm  ee  foil owe t  riret, 
let  r,  be  a  unit  vector  in  the  direction  of  the  epherioel 
ooordinete  system  polsr  axis.  Then  the  magnatio  dipole  field 
X,  -  (Q'p/r1)  <  (3  r  Qos  i  -  r») 

since  r  *  (ix  +  j*y  ♦  fca)/r  where  r  -  (x*  +  y*  ♦  «V  » 

•  •  [  (r*  x  r)  x  r]/sin  e  -  [  (r*rj  r  -  (»*r)  r*]/ein  a,  end  the 
dot  product  r*r#  -  ooe  e,  The  unit  veotor  r(  can  be  expressed 
in  the  rectangular  ooordinstes  by  noting  thst  r,  ■  p/p  end  then 
expressing  p  in  rectangular  ooordinstes.  To  do  this,  let  D 
be  the  depression  angle  of  p  from  the  xy-plane  (the  horiaontal 
plane)  with  0  positive  downward  and  let  a  be  the  direction  of 
p  relative  to  magnetic  north.  Then  the  magnetic  dipole  moment 
p  •  p  (h  cos  n  -  ksin  fl)  in  terms  of  the  unit  vector  k  and 
the  unit  vector  h  -  lain  a  +  j  cos  n  which  has  the  direction 
of  the  horiaontal  component  of  p,  With  those  results,  the  unit 
vector  r0  «  i  (coa  fisin  a)  +  j • (cos  flcos  a)  -  k  sin  n  and 
Hs  -  (c  p/rJ)  ■  [  ( 3/r)  cos  6(ycoa  ♦  -  rain  ♦) 

-  (cos  #’Cos  n  cos  a  sin  S  ein  fl)  ] 
where  cos  e  ■  (1/r)  '(X’cos  n>sin  a  +  ycos  flooe  a  -  a -sin  fl) 
since  cos  e  «  r*r0.  As  can  be  seen  from  this  expression,  for  a 
constant  dipole  moment  magnitude  and  direction  and  a  constant 
earth  field  magnitude  and  direction,  the  magnetic  signal  is  only 
a  function  of  the  rectangular  coordinates  of  the  location  of  the 
magnetometer  relative  to  the  dipole.  In  the  encounter  models, 


22 


I 

i 

|  both  of  thooo  oonditiona  «r«  aatiafiad,  Novavar,  by  allowing  p, 

n,  o,  H,  and  a  to  vary,  tha  axpraaaion  for  H,  ia  applioabla 
to  mora  ganaral  anoountar  modala. 
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Appendix  4.  The  Anderson  Formulation 

In  tho  anoountar  models,  tha  magnetic  aignal  H,  at  a 
sample  point  in  a  atraight  lina  anoountar  oan  ba  rapraaantad  in  a 
fora  daaoribad  by  Anderson  in  Reference  4.  For  oonvanianoa,  tha 
Andaraon  formulation  ia  uaad  in  tha  program  to  datarmina  valuaa 
for  N,.  It  oan  ba  davalopad  aa  follows i  In  Figure  2,  tha 
primad  raotangular  ooordinata  ayatam  ia  auparimpoaad  on  tha 
raotangular  ooordinata  ayatam  of  Figura  1  ao  that  tha  origin  ia 
oolnoidant  with  tha  origin  of  that  ayatam.  A  magnatomatar  ia  in 
a  atraight  lina  anoountar  with  a  magnetic  dipole  that  ia  looatad 
at  tha  origin  of  tha  combined  ayatam.  Tha  combined  ayatam  movea 
with  tha  magnetic  dipola  with  tha  x'-axia  oriented  ao  that  it  ia 
parallel  to  and  in  tha  diroction  p  of  magnatomatar 'a  relative 
to  tha  magnetic  dipola  and  tha  a1 -axis  oriented  ao  that  it  ia 
directed  toward  and  paaaea  through  the  CPA  on  that  track.  Let 
1,  m,  and  n  bo  the  direction  cosines  of  tha  dipole  moment  p 
and  1,,  m,  and  n,  be  the  direction  coainea  of  tha  earth 
magnetic  field  Ht  in  the  primad  ayatam.  Than  tha  unit  vector 
r0  ■*  i '  il  +  j  1  -n  +  A'  -m.  And,  for  pointa  on  tha  relative  track, 
x'  -  a*,  y‘  -  o  and  z '  •  R  where  R  ia  the  alant  range  of  th© 
dipole  at  CPA  and  where  a*  ia  the  algebraic  distance  of  the 
magnetometer  from  CPA  on  the  relative  track.  (It  ia  negative  for 
pointa  before  CPA  and  poaitive  for  pointa  after  CPA.)  Thia 
impliea  that  magnetometer  poaition  vector  in  the  moving 
coordinate  system  is  r  -  i'-fs’/ir)  +  X'fR/r). 
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Magnetic 

East 


Figure  2, 


The  dipole  is  at  the  origin  in  the  unprimed  and  primed 
coordinate  systems.  In  the  primed  coordinate  system, 
the  CPA  is  at  (0,  0,  R)  and,  for  a  time  t  during 
the  encounter,  the  magnetometer  is  at  [s*(t),  0,  R] . 
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From  Appendix  3, 

H,  -  ( c  1  p/r*)  •  ( 3  •  r  •  cob  e  -  r0) . 

With  r  end  r0  expressed  in  terms  of  the  primed  unit  vectors 
end  oos  0  -  r*r0  ■  l- (s’/*)  +  n  •  (R/r) ,  this  becomes 

H,  -  (cp/r*)  •  [  ( 3/r*)  *(l*s»  +  n  *R)  •  (i '  -s'  +  k'  *R) 

-  (i '  •  1  +j '  *m  +  k'  *n) ] . 

Then,  since  H./H,  -  i •  *  1,  +  j'*m,  +  k' *n,  end  H,  -  WH,, 

Ha  ■  (c  *p/r*)  •  [  (2  *1  *  1,  -  m-m,  -  n^n,)*^')' 

+  3  •  (n  •  1,  1  *n,)  •!'  *R  +  (2  *n  *n,  -  1*1,  -  m  *m,)  *R'] . 

The  quantities 

A,-  2  *  1  *  1,  -  m'm,  -  n-n, 

A,  «  3  *  ( n  •  1,  4  1  -n,) 

Aq  -  2  n  -n,  -  1*1,  -  m-m,, 

are  the  Anderson  coefficients.  With  r  -  [(s')*  4  R*]’*  and 

0  «  s'/R#  H,  can  now  be  expressed  as  follows: 

H.(6)  -  (c  p/R3)  T.  hj-F A&) 

where  the  FJU)  -  ft'/(l  +  for  J  -  0,  1,  2  are  the 

Anderson  functions.  This  is  the  Anderson  formulation. 

To  relate  the  Anderson  formulation  for  Hs  to  the 
formulation  for  Hs  in  Appendix  3,  first  note  that 

1  -  r0*i' 
m  -  r0  •  j  ' 
n  -  r0*k' 

and 
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1,  -  (Hj/H,)  -  i* 

®,  "  (VH.)  *j' 
n,  -  (Hg/Hg)  -k'. 

Then  express  r0,  (8,/H,)  and  the  unit  vectors  i',  j '  and  k' 
in  terms  of  the  unit  vectors  i,  j  and  k  and  take  the 
indicated  dot  products.  From  Appendix  3, 

r0  -  i'(cos  ft’Sin  a)  +  j • (cos  n-cos  a)  -  k-sin  n 

and 

hi/hi  ■  j ’cos  *  +  k-sin  ♦  . 

To  express  i*,  j'  and  k'  in  terms  of  the  unit  vectors  i,  j 
and  k,  note  that  the  unprimed  coordinate  system  =an  be 
transformed  to  the  primed  coordinate  system  by  two  rotations  that 
are  defined  as  follows:  First,  rotate  a  coordinate  system  that 
is  coincident  with  the  unprimed  coordinate  system  about  the 
z-axis  through  the  angle  (0  -  it/ 2)  with  positive  angles  clockwise 
(left  hand  rule)  so  that  its  x-axis  is  parallel  to  and  in  the 
direction  of  the  relative  track.  Next,  rotate  this  system  about 
its  x-axis  through  an  angle  6  with  positive  angles 
counterclockwise  (right  hand  rule)  so  that  the  positive  z-axis 
passes  through  the  CPA.  The  angle  6  is  related  to  the 
vertical  separation  z0  of  the  magnetometer  and  the  dipole  and 
the  algebraic  encounter  lateral  range  L  that  is  positive  if  the 
dipole  is  to  the  left  of  the  relative  track.  With  these  sign 
definitions:  L  »  R-sin  6  and  z0  »  R-cos  6.  After  the 
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rotation,  the  auxiliary  coordinate  system  is  coincident  with  the 
primed  coordinate  system. 

These  transformations  can  be  described  in  terms  of  matrix 
equations  as  follows:  Let  (x",yH«zM)  be  the  coordinates  of  a 
point  in  the  coordinate  system  that  is  coincident  with  the 
auxiliary  coordinate  system  after  the  first  rotation.  Then  the 
transformation  from  the  unprimed  coordinate  system  to  this  double 
primed  coordinate  system  is  described  by  the  matrix  equation 

sin  0  cos  0  0  x 

cos  0  sin  00  y 

0  0  1  z 

And  the  transformation  from  the  double  primed  coordinate  system 
to  the  primed  coordinate  system  is  described  by  the  matrix 
equation 

x'  10  0 

y'  *  o  cos  6  sin 

z'  0  -sin  6  cos  6 

Taking  the  product  of  the  rotation  matrices  in  the  indicated 
order  yields  the  matrix  equation 

sin  0  cos  0  Ox 

cos  6-cos  0  cos  6  sin  0  sin  6  y 

sin  S  • cos  0  -sin  6 -sin  <p  cos  6  z 
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y 

z ' 


which  defines  the  transformation  from  the  unprimed  to  the  primed 
coordinate  system.  The  unprimed  vector  components  of  the  unit 
vector  i'  can  be  found  by  transforming  the  coordinates  (1,0,0) 
in  the  primed  system  to  their  corresponding  coordinates  in  the 
unprimed  system  with  the  inverse  of  this  matrix  and  then 
repeating  this  process  for  (0,1,0)  and  (0,0,1)  in  order  to  find 
the  unprimed  unit  vector  components  of  j'  and  Jc'.  However, 
since  the  inverse  transformation  matrix  is  the  transpose  of  this 
matrix,  the  elements  of  the  row  that  corresponds  to  a  primed  unit 
vector  are  the  magnitudes  of  the  unprimed  vectors  that  are  its 
components.  Consequently: 

i '  -  i  -sin  0  +  j  -cos  0 

ji  .  -i -cos -i -cos  0  +  j -cos  i-sin  <p  +  lc-sin  i 
Jc'  -  i-sin  i-cos  0  -  j -sin  6 -sin  0  +  Jc-cos  i. 

Then,  taking  the  dot  products  between  r,  r0  and  these  three 
unit  vectors  as  indicate  above  gives: 

1  =  cos  fl  -cos  (0  -  a) 

m  =  cos  i -cos  n-sin  (0  -  a)  -  sin  i-sin  n 
n  *  -sin  i- cos  n-sin  (0  -  a)  -  cos  i-sin  n 

and 

l,  =  cos  $  -cos  0 

m,  =  cos  i -cos  $-sin  0  -  sin  i-sin  ♦ 
nt  =  -sin  5 -cos  *-sin  0  -  cos  i-sin  #. 

These  are  the  relations  that  are  used  in  the  program  to  determine 
values  for  the  Anderson  coefficients. 
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Appendix  5.  The  Encounter  Equations  of  Motion 

In  the  double  primed  coordinate  system  that  is  defined  in 
Appendix  4,  the  equations  of  motion  of  a  magnetometer  relative  to 
a  submarine  (dipole)  are: 

x"(t)  -  s' (t) 
y"(t)  -  -l 

Z"(t)  -  20 

where  L  is  the  algebraic  encounter  lateral  range  that  is 
defined  in  Appendix  4,  z0  is  the  vertical  separation  between  the 
magnetometer  and  the  submarine  and  s'(t)  is  the  distance  of  the 
magnetometer  from  the  CPA  on  the  relative  track.  With  w  the 
speed  of  the  magnetometer  relative  to  the  submarine  and  t  a 
relative  time  parameter,  s'(t)  «=  wt  .  These  equations  can  be 
considered  to  be  the  ones  used  in  the  program  to  describe  the 
motion  of  a  magnetometer  relative  to  a  submarine.  There,  t  is 
determined  by  t  =  [J  -  (M-l)/2]’6t  where  the  index  J  *  1,  2, 
*•*,  M  and  «5t,  a  time  step,  is  the  time  between  samples.  Note, 
when  t  =  0,  the  magnetometer  is  at  the  CPA. 

In  the  coordinate  system  of  Figure  1  in  a  straight  line 
encounter  as  defined  in  the  encounter  models,  the  equations  of 
motion  of  a  magnetometer  relative  to  a  submarine  can  be  written 
as  follows: 


x(t) 

-  s'(t) 

•sin  <p  + 

L'COS  0 

y(t) 

-  s'(t) 

•cos  0 

L -sin  0 

z(t) 

“ 
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since  the  transformation  from  the  double  primed  coordinate  system 
to  the  primed  coordinate  system  is  determined  by  the  matrix 
equation 


X 

sin  0  -cos  0  0 

X" 

y 

= 

cos  0  sin  0  0 

y " 

z 

0  0  1 

z» 

x"  *  s"(t)  and  y"  ■  -L.  The  above  equations  and  the  expression 
for  Hs  in  Appendix  3  could  have  been  used  in  the  program  to 
evaluate  the  magnetic  signal.  In  particular,  with  these 
equations  of  motion  and  with  the  two  relations  L  =>  R-sin  5  and 
z0  =  R'cos  6,  the  expression  for  Hs  in  Appendix  3  can  be  written 
in  terms  of  n,  a,  t,  <p,  6  and  R  so  that  it  is  identical  in 
appearance  to  the  Anderson  formulation  for  Hs  in  terms  of  these 
quantities.  The  definition  of  <5  in  Appendix  4  in  terms  of  a 
counterclockwise  rotation  results  in  a  definition  of  the 
algebraic  lateral  range  that  is  consistent  with  some  that  have 
been  used  elsewhere. 

In  the  program,  the  relation  w  =  v  -  u  is  the  basis  for 
determining  the  relative  speed  w.  In  this  relation,  v  is  the 
velocity  of  the  magnetometer,  u  is  the  velocity  of  the  submarine 
(dipole)  and  w  is  the  velocity  of  the  magnetometer  relative  to 
the  submarine.  This  relation  implies  the  following  equations: 
w„  -  V,  -  u„,  wy  -  vy  -  uy  and  w,  =  v,  -  u,  where  the  coordinates 
x,  y  and  z  refer  to  a  fixed  coordinate  system  with  the  same 
orientation  as  that  of  Figure  1.  In  the  encounters  of  the 
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models,  v,  -  v*sin  a,  vy  -  v*cos  a  and  v,  ■  0  where  a  is  the 
magnetic  course  and  v  is  speed  of  the  magnetometer.  And,  in 
addition,  u,  -  u-sin  R,  Uy  -  u*cos  R  and  u,  -  0  where  R  is 
the  magnetic  course  and  u  is  the  speed  of  the  submarine.  The 
relative  magnetic  course  <f>  and  the  relative  speed  of  the 
magnetometer  are  defined  by  w,  ■  wsin  <t>  and  wy  -  w*cos  <fi.  In 
the  program,  <p  and  w  are  determined  with  a  rectangular  to 
polar  conversion  routine  where  w  -  (  w„a  +  wy*),/8  and  where  <p 
is  determined  by  sin'1  (w*/w)  and  cos'(wy/w). 
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Appendix  6.  Th*  Submaria*  Magnetic  Dipol* 

In  the  encounter  models,  the  magnitude  and  direction  of  a 
submarine's  dipole  moment  p  are  determined  by  first  determining 
its  components  in  the  rectangular  coordinate  system  of  Figure  1. 
In  that  coordinate  system, 

P«  -  (Plp  +  Pu)  -sin  fi  +  (Ptf  +  Pn)  *cos  6 

Py  -  (Plp  +  Pu)  *cos  6  -  (Ptp  +  Pn)  ' sin  fl 

P,  -  "(Pvp  +  Pv) 

where  fi  is  the  submarine's  magnetic  course,  p^,  p^  and  p^ 
are  the  permanent  and  pu,  pfl  and  pw  are  the  induced 
longitudinal,  transverse  and  vertical  magnetic  dipole  moments  of 
the  submarine.  These  relations  are  based  on  the  following  sign 
convention: 

pL  is  positive  when  pL  is  directed  from  stern  to  bow. 

pT  is  positive  when  pT  is  directed  from  port  to  starboard. 

pv  is  positive  when  pv  is  directed  downward. 

The  permanent  dipole  moments  are  input  parameters  in  the  program 
and  the  induced  dipole  moments  are  determined  in  a  way  that  is 
similar  to  one  that  is  described  in  Reference  1.  In  the 
encounter  models,  a  submarine  is  a  ferromagnetic  prolate 
ellipsoid  with  the  major  axis  the  submarine's  longitudinal  axis 
and  the  equal  minor  axes  the  submarine's  transverse  and  vertical 
axes.  And,  the  induced  dipole  moments  are: 

Pu  =  ’^EL 
Pn  =  ’®et 
Pv,  =  K  ' 
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where  H^,  and  are  the  vector  components  of  the  earth 

magnetic  field  in  the  rectangular  coordinate  system  defined  by 
the  ellipsoid  axes  and  the  submarine  magnetic  moment  sign 
convention.  The  magnitude  of  these  vector  components  are: 

Hel  =  He-cos  * -cos  B 
Her  =  -HE-cos  *-sin  6 
Hp,  =  HE  -  sin  *. 

The  earth  magnetic  field  dip  angle  *  and  the  earth  magnetic 
field  magnitude  HE  can  each  be  chosen  to  be  input  parameters  in 
the  program  or  the  can  be  computed  by  the  program  as  described  in 
Appendix  7.  (Values  for  *  and  HE  can  be  found  from  magnetic 
data  charts,  for  example,  see  References  5  and  References  6.) 

By  using  the  above  relations, 

p,  =  He*cos  #-(kL  -  kT)  <sin  B-cos  B  +  (p^-sin  B  +  p^-cos  B) 
py  =  HE*cos  #-(kL-cosJB  +  kT -sin'S)  +  (p^-cos  B  -  p^-sin  B) 

P,  =  -(iWsin  *  +  Pyp) 

In  the  encounter  models,  the  values  of  the  permeability 
coefficients  kT  and  ky  are  related  to  submarine 

displacement  by  the  following  relations: 
kL  *  fL*W 
kT  =  fT-W 
K  m  fv-w 

where  W  is  the  submarine  displacement  in  tons  and  fL,  fT  and 
fv  are  permeability  factors  that  are  determined  by  a  submarine's 
magnetic  characteristics.  In  the  program,  the  units  of  HE  are 
gamma,  the  units  of  p  are  in  oersted-centimeter3,  the  units  of 
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k  ara  oersted-centimetar'/gamma  and  tha  units  of  f  ara 
oersted-cant  imeter'/gamma-ton.  If  tha  units  of  Hj  wars  gamma , 
but  tha  units  of  p  wars  gamma-foot',  than  tha  units  of  k 
would  ba  foot'  and  tha  units  of  f  would  ba  foot'/ ton.  To 
convart  p  in  gamma- foot'  to  oarstad-oantimatar*  or  k  in  foot* 
to  oarstad-cantimatar'/gamma  or  f  in  foot'/ton  to  oaratad- 
cantimatar'/gamma-ton,  divida  by  3. S3.  (Tha  program  dafault 
values  for  fL,  fT  and  fv  ara  valuas  from  Rafaranoa  1  in 
foot'/ton  that  have  been  divided  by  3.53  to  give  valuas  in 
oersted-centimeter’/gamma-ton. )  To  oonvart  p  in  webar-metar  to 
oersted-centimeter’,  multiply  by  [ 1/ (4  •  w ) ]  •  10"  and  to  convart 
p  in  ampere-meter*  to  oersted-centimeter',  multiply  by  1*10'. 
Note,  He  in  oersted  equals  10 8  times  Ha  in  gamma,  H,  in  gamma 
equals  HE  in  nanotesla  and  HE  in  tesla  equals  Hc  in 
weber/meter',  so  HE  in  gamma  equals  10*  times  Ht  in 
weber/meter*  and  HE  in  oersted  equals  104  times  HE  in 
weber/meter*. 


Appendix  7.  The  Perth  Magnetic  Pleld  mi  Dtp  Angie 

An  Auxiliary  magnetic  field  model  it  described  in  thia 
appendix.  Tha  modal  ia  tha  baaia  for  a  dafault  ohoiea  for  aithar 
tha  valua  of  tha  aarth  magnetic  fiald  magnitude  parameter  H,  or 
tha  aarth  magnatio  fiald  dip  angla  paramatar  ♦.  Relative  to 
anoountar  modal  accuracy,  tha  dafault  valuaa  ahould  ba  adaquata 
in  moat  oaaaa. 

In  tha  modal ,  tha  a^rth  magnatio  fiald  ia  ganaratad  by  a 
magnatio  dipola  that  ia  looatad  at  tha  earth's  oantar  and  tha 
aarth  ia  a  nonmagnatio  aphara  of  radiua  r*.  Nith  p*  tha 
magiutuda  of  tha  dipola  moment,  tha  magnituda  of  tha  aarth  fiald 
at  a  point  ia 

Ht  -  (Hf  *Mi) '  *  -  (o*p,/rJ)  •  (3  coa'e  ♦  l)'* 
whara  G  ia  tha  polar  angla  of  tha  point  in  a  apharical 
coordinate  system  and  c  ia  a  conatant  whose  valua  ia  determined 
by  the  choice  of  unite.  The  dipola  moment  ia  coincident  with  and 
in  the  direction  of  the  polar  axia  which  ia  directed  toward  the 
earth'a  aouthern  hemisphere.  In  thia  coordinate  system,  at  any 
point  on  the  surface  of  the  earth: 

H{  -  H,0  ■  (3  >008*0  +  I)’  * 

where  Ht0  ia  the  value  of  Hf  at  the  magnetic  equator  which  is 
defined  by  the  points  on  the  earth  where  0  -  90*.  In  terms  of 
the  dip  angle,  at  any  point  on  the  earth's  surface: 

HE  -  2  -Ht0  •  (3  cos**  +  1)  ’ 
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This  expression  ean  tea  obtained  fey  noting  that  #  oan  ba 
defined  in  terms  of  tha  or  tha  r  and  t  oomponanta  of  N,  in 
tha  apharical  uoordinataa  aa  follows: 

ain  ♦  -  -H^N,  -  -2  'H*,'  (ooa  •)/!«« 
and 


cos  ♦  »  *WHt  -  Hw’  (ain  •)/«,. 

Based  on  thaaa  relations, 

(  (ain  ♦ ) / a ] #  i  ooa'#  -  (HM)‘ 
and 

H,  *  Htu  ■  ( (ain*)/4  >  Goa'#]'*  •  2  Hjjo  (3  coa*a  4  l)'* 

Tha  dip  angle  iw  determined  from  a  magnetic  latitude  for  tha 
encounter  region, 

Tha  magnetic  latitude  and  longitude  of  the  representative 
point  of  an  encounter  region  can  ba  defined  in  terms  of  its 
geographic  latitude  and  longitude  by  using  the  following 
transformations:  First,  convert  the  latitude  and  longitude  of 

the  point  to  rectangular  coordinates  in  a  right-handed  coordinate 
system  whose  origin  is  at  the  center  of  a  spherical  earth,  whose 
positive  E-axis  passes  through  its  north  geographic  pole  and 
whose  positive  y-axis  passes  through  the  point  on  its  equator 
with  latitude  0*  and  longitude  O’.  Next,  rotate  a  coordinate 
system  that  is  coincident  with  this  system  in  a  clockwise 
direction  (left  hand  rule)  about  its  z-axis  so  that  its  positive 
y-axis  passes  through  the  point  with  latitude  O’  and  with 
longitude  equal  to  that  of  the  north  magnetic  pole.  Then,  rotate 
the  system  in  a  clockwise  direction  about  its  x-axis  so  that  Its 
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positive  i-axia  paaaaa  through  tha  north  msgnatio  polo.  Next, 
transform  tha  raotangular  ooordinatos  of  tha  rapraaantativa  point 
in  thia  system  to  ita  ooordinataa  in  tha  apharieal  ooordinata 
ayatan  that  ia  aaaoeiatad  with  it.  Than,  with  n  tha  polar 
angla  of  tha  rapraaantativa  point  in  thia  ayatan,  tha  msgnatio 
latituda  of  tha  rapraaantativa  point  ia  X*  -  SO*  -  n.  tinea 
tan  •  -  -loot  a  from  abova  and  •  -  180 ’  -  n,  tha  dip  angla  # 
ia  givan  by  tha  following  ralationi 
•  ■  tan '(2  -  tan  L*) . 

In  tha  program,  tha  tranaformationa  daacribad  abova  ara 
accompliahad  in  part  by  a  raotangular  to  polar  oonvaraion 
routina.  In  particular,  by  rotating  tha  final  raotangular 
ooordinata  aystam  about  ita  z-axia  so  chat  tha  x-coordinata  ia 
aaro,  tha  polar  angla  of  tha  polar  ooordinataa  of  tha 
raprasentative  point  in  tha  raaulting  yz-plana  datarminaa  I*. 

Soma  values  of  *  and  Ht  that  are  listed  in  Table  1  were 
generated  using  tha  program.  Tha  latituda  and  longitude  of  the 
magnetic  pole  that  are  in  the  program  are  76*  N  and  100*  w. 
Values  of  *  from  Reference  6  and  of  HE  from  Peference  7  are 
also  listed  in  Table  1.  A  comparison  of  the  program  values  with 
these  values  gives  an  indication  of  the  errors  inherent  in  the 
procedure. 
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Encounter 

Encounter 

Program 

Chart  130 

Program 

Chart 

Latitude 

Longitude 

Dip  6 

Dip  • 

Field  H« 

Field 

60  *  N 

180  ’  W 

74* 

70* 

.63 

.83 

60'N 

180*W 

78* 

73* 

.66 

.85 

«0*N 

90’W 

83’ 

84* 

.68 

.61 

60*N 

30*W 

78* 

74* 

.64 

.52 

30*N 

150*W 

88* 

SO* 

.81 

.41 

30’N 

90  *W 

62  * 

60* 

.88 

.51 

30*N 

60*W 

89* 

89* 

.83 

.48 

30  'N 

30*E 

37* 

43* 

.41 

.42 

30  *  N 

60  *  E 

31* 

45* 

.39 

.46 

30*N 

90  ‘  E 

30* 

44* 

.39 

.44 

30  *N 

150*E 

42* 

40* 

.43 

.40 

0* 

90  *W 

26* 

18* 

.38 

.32 

O’ 

30  *W 

9* 

-5* 

.35 

.28 

0* 

0* 

-48* 

-25* 

.35 

.31 

0* 

60  *  E 

-25* 

-20* 

.38 

.37 

30  *S 

180  *W 

-45* 

-55* 

.44 

.50 

30  ’  S 

90 "  W 

-30* 

-34* 

.  39 

.31 

30  ‘ S 

90 '  E 

-62  ' 

-66* 

.55 

.52 

60 ’S 

3  0  *  W 

-69’ 

-56* 

.60 

.34 

m 

H, 


Table  1.  Some  corresponding  values  for  the  encounter  dip  angle 
*  and  the  earth  magnetic  field  intensity  Hf  in 
oersted.  Chart  #33  refers  to  Reference  5  and  Chart  #39 
refers  to  Reference  6. 


A  small  difference  between  a  chart  value  and  a  computed  value  for 
the  dip  angle  *  at  some  location  does  not  imply  a  small 
difference  between  the  chart  value  and  the  computed  value  for  the 
earth  magnetic  field  intensity  at  that  location.  For  example, 
note  the  values  for  60 ‘N  latitude  and  30 'W  longitude. 
(Positive  values  of  $  indicate  the  inclination  or  dip  below  the 
horizontal  of  the  north  seeking  end  of  a  dip  needle.  Negative 
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valuta  indicate  tha  inclination  below  the  horizontal  of  tha  aouth 
■••king  and.) 

An  axtanaion  of  tha  above  prooadura  oan  be  made  for  finding 
tha  nagnetio  variation  at  a  location.  However,  tha  valuaa 
generated  by  uaing  tha  prooadura  are  generally  unaatiafaotory. 
Nagnatio  variation  valuaa  are  ohartad  in  Reference  7, 


.>pp«ndtx  «.  An  Alternative  inoounter  Model 

The  alternative  model  that  is  rsfsrrsd  to  in  Ssotion  HZ  la 
described  in  mors  detail  in  this  appendix.  In  ths  nodal,  ths 
dataotion  range  is  tha  slant  ranga  R  at  ths  CPA  for  an 
anoountar  with  a  spaoifiad  dataotion  probability  (usually  .5)  is 
dafinad  by: 

R  -  [C  *p/H,]  ui 

whara  c  and  p  are  dofinad  in  Appandix  3  and  H,  raprasants 
a  minimum  detectable  average  nagnetic  signal  that  is  defined  by: 

Hs  -  (ORF) -Nm 

where  ORF  is  a  signal-to-noise  ratio  called  the  operator 
recognition  factor  and  Nw  represents  the  nagnetic  noise. 
Combining  these  two  relations  gives: 

R  -  (c-p/[  (ORF)  -NM]  )u\ 

The  value  for  ORF  depends  both  on  the  specified  detection 
probability  and  on  a  specified  or  implied  false  alarm 
probability. 

The  Anderson  formulation  is  consistent  with  these  relations 
in  an  approximate  sense  if  the  average  magnetic  signal  H  is 
defined  as  a  root  mean  square  value  such  that  H  =  (c-p/R3)  -k 
where  k  is  an  encounter  parameter  defined  by: 
k  =  {Z  [Z  A,  •  Fj ( Bk)  ]?)’/2 

with  the  first  sum  index  K  =  1,  2,  •••  ,  M  and  the  second  sum 

index  J  *  0,  1,  2.  For  a  particular  encounter  geometry,  k  is 
constant  and  this  suggests  that  the  two  encounter  models  could  be 
used  to  determine  an  average  value  for  k  for  an  encounter 
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region  baaad  on  avaraga  aubmarina  magnatic  oharactariatioa. 

Valuaa  for  both  k  and  H  ara  ganaratad  by  tha  program  and  ouch 
valuaa  can  giva  an  indioation  of  tha  magnituda  of  tha  diffaranoaa 
in  dataotion  ranga  aatimataa  that  ara  baaad  on  thia  modal  and 
thoaa  that  ara  baaad  on  aithar  of  tha  othar  two  anoountar  modal a. 
A  mora  datailad  compariaon  of  thaaa  ancountar  modala  ia  daaoribad 
in  Rafaranoa  8. 
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Appendix  I.  An  sxnaple  of  the  Program  output 

Tho  program  in  dssigntd  to  gonorato  valuta  for  tha  following 
quantities!  ancountar  paramatsrs,  lataral  range  funotiona  for  tha 
oroaaoorrelation  encounter  model  and  for  the  aquare  law  encounter 
model,  average  magnatio  aignal,  alant  range  at  CPA,  encounter 
parameter  valuta  and  magnetic  aignal  valuta.  Tha  valuta  can  be 
aaved  as  a  program  data  file  and/or  they  can  be  printed.  Tables 
2  through  7  are  examples  of  the  program's  printed  output. 

Figures  3  and  5  are  plots  of  tha  lateral  range  function 
values  that  are  listed  in  Tables  3  and  6  for  the  square  law 
encounter  model.  For  either  the  square  law  or  the 
crosscorrelation  encounter  model,  as  the  magnitude  of  the  lateral 
range  increases,  lateral  range  function  values  (detection 
probabilities)  do  not  approach  zero  but  approach  the  value  of  the 
encounter  false  alarm  probability.  A  necessary  condition  for  the 
formal  definition  of  sweep  width  given  in  Reference  9  to  be 
meaningful  is  that  the  lateral  range  function  approach  zero  as 
the  lateral  range  increases  without  limit.  Since  this  is  not  the 
case  for  either  encounter  model,  the  formal  definition  of  sweep 
width  must  be  modified  with  their  use.  Reference  10  provides  an 
example  of  this. 

Figures  4  and  6  are  plots  of  the  magnetic  signal  values  that 
are  listed  in  Tables  4  and  7. 
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program  fila  nams 
program  data  file  nama 
magnatio  data  til  a  nama 
prooaaaing  data  fila  nama 
Kinamatic  data  fila  nama 

oombinad  magnatio,  processing  &  kinamatic  data  film  nama 

ancountar  latituda  (daeimal  dagraaa) 

anoountar  lonoituda  (daoimal  dagraaa) 

anoountar  variation  (daoimal  dagraaa) 

anoountar  dip  anola  (daoimal  daoraaa) 

anoountar  magnatio  fiald  intanaity  (gamma) 

parmanant  longitudinal  moraant  (oaratad-oml) 

parmanant  tranavarsa  momant  (oaratad-c*J) 

parmanant  vartioal  momant  (oaratad-omS) 

targat  diaplaoamant  (tona) 

targat  longitudinal  parmaability  ooaffioiant 

targat  tranavaraa  parmaability  ooaffioiant 

targat  vartioal  parmaability  ooaffioiant 

targat  longitudinal  parmaability  factor 

targat  tranavaraa  parmaability  factor 

targat  vartioal  parmaability  factor 

sampling  pariod  (saconda) 

intagration  tima  (saconda) 

adjustad  intagration  tima  (saconda) 

numbar  of  samplas  par  ancountar 

magnatomatar  coursa  (daeimal  dagraaa) 

magnatomstar  spaad  (knots) 

magnatomatar  altituda  (matars) 

targat  coursa  (decimal  dagraas) 

targat  spead  (knots) 

targat  depth  (maters) 

magnetometer  relative  course  (decimal  degrees) 

magnetometer  relative  spaed  (knots) 

magnetometer-target  vertical  separation  (matars) 

target  induced  longitudinal  dipole  momant  (oersted-cm3) 

target  induced  transverse  dipole  momant  (oersted-cm3 ) 

target  induced  vertical  dipole  moment  (oersted-cmi ) 

magnetic  dipole  moment  (oersted-cm3) 

dipole  moment  azimuth  (decimal  degrees) 

dipole  moment  depression  angle  (decimal  degrees) 

distance  between  samples  on  the  relative  track  (meters) 

false  alarm  rate  (false  alarms  per  hour) 

false  alarm  probability 

magnetic  noise  (gamma) 

maximum  lateral  range  (meters) 

lateral  range  step  (meters) 

number  of  lateral  range  function  values 


HAD. BAS 
data. mad 
data. mag 
data. pro 
data. kin 
data.mpk 
49 
-45 
-85 

68.41608 

58035.05 

0 

0 

0 

1800 

13140 

3880 

3880 

7.3 

1.6 

1.6 

.5 

30 

30.5 

41 

290 

180 

100 

20 

10 

100 

286.8202 

180.2776 

200 

2 . 017795E+08 
-4.422564E+07 
1 . 58  099E+08 
2 . 601272E+08 
7.637512 
37.42885 
46.3714 
3 

1 . 666667E-02 
.35 
500 
20 
51 


Table  2. 


An  example  of  an  encounter  parameter  values  printout 


,1.1'.  1 . '  1- 

1  l.i*.  t-i.ii  i.in.i<3 

!  Ufi.vtii.fi  value 

ti 

L 

■) 

p(»l) 

H 

R 

k 

meters 

gamma 

matara 

•■500 

.  2  1559 6$ 

. 0255824 

.  184247 

538 . 5165 

2.818486 

-<180 

.2721803 

. 029692 

.21034 

520 

2.879916 

-460 

.3464673 

3 . 582633E-02 

.2410144 

501.5974 

2,943744 

-4  4  0 

.4415205 

4 . 52  9992E-02 

.2785504 

483.3218 

3.009938 

-4  20 

.  5573696 

6 , 047795E-02 

.3231911 

465.1681 

3 . 078418 

-400 

. 6872897 

.0856887 

.3764825 

447.2136 

3 . 149025 

-380 

.8146549 

. 1287804 

.4404623 

429.4182 

3.221497 

-360 

.9159235 

.2030898 

.5177304 

411.8252 

3.29544 

-340 

.9745625 

. 3272475 

.6108596 

394.4617 

3 .370274 

-320 

.9958122 

. 5141411 

.7233319 

377.3593 

3.445185 

-300 

. 9997211 

.7388175 

.8597726 

3C0.5551 

3 . 519032 

-2  8  0 

. 9999951 

.9177597 

1.025158 

344 . 093 

3 . 590261 

-260 

1 

.9896971 

1.224529 

328.0244 

3 . 656775 

-240 

1 

. 999714 

1.463644 

312.41 

3 . 71578 

-220 

1 

. 9999992 

1.74981 

297.3214 

3 . 763606 

-200 

1 

1 

2.065394 

282.8427 

3 . 795517 

-180 

1 

1 

2.470496 

269.0725 

3 . 805534 

-160 

1 

1 

2 , 901352 

256.125 

3.786234 

-140 

1 

1 

3.359555 

244 . 1311 

3 . 72898 

-120 

1 

1 

3.819509 

233.2361 

3 . 624058 

-100 

1 

1 

4.239317 

223.6068 

3 . 461596 

-80 

1 

1 

4 .532758 

215.4066 

3 . 232958 

-60 

1 

1 

4.618943 

2CB.8061 

2 . 932828 

-40 

1 

1 

4.44342 

203 .9608 

2.5617 

-20 

1 

1 

3.97739 

200.9975 

2.128285 

0 

1 

1 

3.332139 

200 

1 . 651697 

2  0 

1 

,  9  Ci  0  Q  Q  9  3 

2.540666 

200.9975 

1.165581 

40 

0  q  n  r.  -  • 

.  9  06  7 1 3 

1 .76834 

203 .9608 

.7394173 

60 

.9933111 

.467368 

1 .110188 

208 .8061 

.5637167 

0  0 

.  9  9  9  S  7  j  3 

.7874625 

.9039734 

215.4066 

.7782836 

3  0, 

,9999009 

.9856169 

1 .086971 

223 . 6068 

1.134662 

120 

1 

.9992399 

1.254863 

233.2381 

1 . 490748 

140 

1 

.9998478 

1.308583 

244 . 1311 

1.811734 

160 

1 

. 9998431 

1.29482 

256.125 

2 . 089946 

180 

1 

.9994955 

1.230946 

269.0725 

2 . 326072 

200 

1 

.9974561 

1.131703 

282.8427 

2 . 523877 

22  3 

1 

.9872744 

1.017684 

297 . 3214 

2 . 688154 

24  0 

.999999 

.9510013 

.9018685 

312.41 

2.823785 

260 

.  9  9  9 f ■  7  4  7 

.8646194 

.  7917129 

328 . 0244 

2.9353C7 

2  £  7- 

.9996544 

.7242736 

.6909791 

344 . 093 

3 . 026728 

3  e  c 

.5973502 

.5582789 

.6011183 

360.5551 

3 . 10149 

320 

.9876337 

.4037092 

.5222269 

377 .3593 

3 . 162495 

34  0 

.9C08865 

.2816961 

.4536729 

394.4617 

3.212156 

3  ( 

. 9086 40" 

.1948741 

.  394481  1 

411.8252 

3 . 252463 

3  F  i 

.8301703 

.1365131 

.3443362 

429.4182 

3 . 285053 

4  ■: : 

. 7333005 

9 . 8204  3  IE-02 

.3014853 

447.2136 

3.311281 

4  2  3 

.  6  2  9  9  3  8  7 

. 0731037 

.2645045 

465. 1881 

3.332228 

4  4  . 

.  5  3  0  2  "  0  4 

5.6466368-02 

.2325951 

483.3218 

3. 348799 

4  6  0 

.440881“ 

4 . 5229331  -02 

.205043 

501 . 5974 

3.361725 

4?.) 

.  j  6  4  !■)  4  4  f ' 

3 . 7471O6E-02 

.1812233 

520 

3. 373603 

50  3 

.  3013791 

3  .  199235F-02 

. 1605954 

538. 5165 

3 . 378925 

Table  3.  An  example  of  a  lateral  range  function  values  printout. 
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find 


heading  for  the 
the  heading  for 


crosscorrelation  values  is  p(cc) 
the  square  law  values  is  p(sl). 


45 


Lateral  ^ange 


Figure  3.  A  plot  of  the  square  law  lateral  range  function  values 
that  are  listed  in  Table  3.  The  horizontal  axis  is 
encounter  horizontal  range  at  CPA  in  meters.  The 
vertical  axis  is  encounter  detection  probability. 


aft 


data. mad  magnetic  signal  values 

for  a  lateral  range  of  40  meters 

magnetic  signal 
Hs  in  gamma 


relative  CPA  distance 
s  in  meters 


-927.4:62 
-881.0568 
-834 . 6854 
-788.314 
-741.9426 
-695. 5712 
-649.1998 
-602.8283 
-556.4569 
-510.0855 
-463.7141 
-417.3427 
-370.9713 
-324 . 5999 
-278.2285 
-231.8571 
-185.4856 
-139.1142 
-92.74282 
-46.37141 
0 

46.37141 

92 . 74282 

139. 1 142 

185.4856 

231.8571 

278.2285 

324 . 599? 

370.971? 

417.3427 

463.7141 

510.0855 

556.456° 

602.8283 

649.1998 

695.5712 

741.9426 

788.314 

834.6854 

661 .0568 

927.4262 


-1 . 619 177E-02 
-1.8343 1 5E-02 
-2 . 086144E-02 
-2.38 1637E-02 
-2 . 728804E-02 
-.0313641 
-3 . 613118E-02 
-4 . 1654 12E-02 
-4 . 792992E-02 
— 5.478963E-02 
-6. 169327E-02 
-6 . 730653E-02 
-.0686387 
-5. 933568E-02 
-2. 652314E-02 
.0540608 
.2215785 
.5163918 
.9104893 
1.184566 
1 . 003207 
. 3837067 
-.2252113 
-. 5323231 
-.5837545 
- . 5168925 
-.4202569 
-  .3304064 
-.2572638 
-.2007385 
-.1577988 
-.125285 
-.1005532 
-8. 158605E-02 
-6. 689598E-02 
-5. 539922E-02 
-4.63  07 19E-02 
-3 . 9043625-02 
-5 . 3 18428L-02 
-.0284142 
-2 . 449729E-02 


Table  4.  An  example  of  a  magnetic  signal  values  printout. 
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I 


Figure  4.  A  plot  of  the  magnetic  signal  values  listed  in  Table  4. 

The  horizontal  axis  is  relative  CPA  distance  in  meters. 
The  vertical  axis  is  magnetic  signal  in  gamma. 
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program  fila  name 
program  data  til*  nama 
magnetic  data  file  nama 
processing  data  file  nama 
kinematic  data  fila  nama 

combined  magnatic.  processing  I  kinamatio  data  file  nama 

ancountar  iatituda  (daoimal  degrees) 

ancountar  longitude  (daoimal  dagraaa) 

ancountar  variation  (daoimal  dagraaa) 

ancountar  dip  angla  (daoimal  dagraaa) 

anoountar  magnetic  fiald  intanaity  (gamma) 

permanent  longitudinal  moment  (oaratad-oml ) 

permanent  tranavaraa  moment  (oaratad-ca3) 

permanent  vartioal  moment  (oerated-aal) 

target  diaplaeamant  (tone) 

target  longitudinal  permeability  ooaffioient 

target  tranavaraa  permeability  ooaffioient 

target  vertical  permeability  ooaffioient 

target  longitudinal  permeability  faotor 

target  tranavaraa  permeability  faotor 

target  vertical  permeability  factor 

aampling  period  (saeonds) 

integration  time  (aaoonda) 

adjuated  integration  time  (aaoonda) 

number  of  aaoplea  par  anoountar 

magnetometer  couraa  (decimal  dagraaa) 

magnetometer  apaad  (knota) 

magnetometer  altitude  (metara) 

target  couraa  (decimal  dagraaa) 

target  spaed  (knota) 

target  depth  (metara) 

magnetometer  relative  couraa  (deoiaal  dagraaa) 

magnetometer  relative  speed  (knots) 

magnetometer-target  vertical  separation  (maters) 

target  induced  longitudinal  dipole  moment  (oersted-ami ) 

target  induced  tranavaraa  dipole  moment  ( oersted-, ?ml ) 

target  induced  vertical  dipola  moment  (oaratad-cml) 

magnetic  dipola  moment  (oaratad-cml) 

dipola  moment  azimuth  (decimal  dagraaa) 

dipola  moment  depression  angla  (decimal  degrees) 

distance  between  samples  on  the  relative  track  (meters) 

false  alarm  rata  (falaa  alarms  par  hour) 

false  alarm  probability 

magnetic  noise  (gamma) 

maximum  lateral  range  (meters) 

lateral  range  step  (matara) 

number  of  lateral  range  function  values 


NAO . hAd 
data i . ma.l 

data. meg 
data, pro 
detal.kin 
datal.mpk 
41 
•U 
-3» 

••.41*01 

19018,08 

0 

0 

0 

1100 

13140 

3M0 

aaao 

7.3 

i.a 

l.e 

,» 

30 

30.  S 
41 

340 

lao 

100 

373 

10 

100 

390.9931 

l'Q.462 

200 

i.339*aas*oa 
3. S22131E*07 
l,»«099r*08 
a  .  14  4  5*4  E  +  08 
295.4033 
47,49386 
41.84*62 
1 

1.S66667E-02 

.  13 
300 

20 

3 1 


Table  5.  An  example  of  an  encounter  parameter  values  printout. 
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UMMl  MM*  UMloduM  \M\i96 


u 

9(*\) 

H 

4 

4 

Mt*r* 

44M4 

i*i«n 

9,1  IMMS'Oi 

1, 41140*l-0| 

,419441) 

9)4,1141 

4.41*14* 

-4*0 

.  1041019 

4,441 t*4K*0i 

, 4201412 

120 

4,0*419* 

-490 

. 14481*1 

1  .97*4048*01 

,441*441 

904,4414 

4 , 0 1  *  V*  7 

*44« 

,  1449141 

1 , 094IO4R-0I 

, 441,14  2 

491,  m* 

\,»71t4) 

*42>' 

.  1101*11 

2 , 2644k  78-01 

,  1911904 

4*4.4991 

1 , 4)0) 

-400 

,2011041 

1 , 4110044*01 

,1444494 

441,1414 

1,49414 

-)90 

.1441191 

1,1192941-01 

,1442411 

414,4491 

1,14491 2 

•)*0 

.2012422 

2 ,3244148-01 

,2914313 

444,1211 

1,44149* 

-540 

,2444001 

1 , 9111041-01 

,2444034 

394,4441 

1,42*114 

-  J*0 

,4444244 

1.0414191-01 

,4041014 

311,219) 

1,4409** 

•)00 

•444449 

1,0941411-01 

,4430444 

140, ill) 

1,41901* 

-.'10 

.  1*442044 

. 4092111 

.1192143 

144,092 

1.19409 

-240 

,4944444 

,1910014 

,4943413 

314,0144 

1,9411*9 

•240 

,9114141 

,1114411 

,4440929 

111,44 

4.01904) 

*220 

,44^2*14 

,1944401 

1,031113 

191,1144 

4.109*9) 

-100 

,9994411 

,4339119 

4, 141044 

141.9411 

4.311441 

-190 

,9949999 

,9144144 

4,141114 

149,0111 

4,341429 

•140 

1 

,9494997 

4,904141 

134,114 

4,41*091 

-140 

1 

,9949949 

1, 341414 

144.4244 

1,4)4011 

•110 

V 

4 

1,174494 

133,1344 

1,141399 

•100 

1 

4 

1,144494 

113,4044 

1,414014 

-«Q 

1 

1 

1,944111 

119,4044 

1,419009 

•  40 

1 

1 

4, 144407 

104.9091 

1,040**9 

•40 

1 

t 

4 , 7 14QI4 

10) , 9*09 

3,0941*4 

-20 

1 

l 

4,9214)2 

200,9979 

3,009993 

0 

V 

l 

4**3249 

300 

«,*04’*9 

20 

l 

l 

4,944144 

200,9*79 

2.742*97 

40 

4 

4 

4,04414 

901.9*04 

4, *34*04 

40 

1 

1 

1,4*1993 

209,90*1 

3 . 29(933 

10 

1 

1 

2.414174 

119,40** 

9.04926) 

100 

1 

, 9*994*1 

2. 141797 

22). *0*9 

1,91221) 

120 

1 

,4*09414 

1,441014 

111.2)91 

1 , *09*3 

140 

. 9994412 

.  4493344 

1,4423*2 

144.1311 

1,44**91 

UO 

.9911449 

.94791 

4.440499 

19*. 423 

1,349197 

110 

,9773194 

.1194244 

.9491192 

3*9,0723 

1.31299 

200 

,9211947 

. 2212309 

,7337172 

242.4427 

1.334)14 

220 

,••'14791 

,1414171 

.9979109 

2*7.3314 

1,39*3* 

240 

,9071171 

,  13*1274 

.4944332 

711.41 

1,4*099 

240 

.7414247 

,  100*447 

.4141041 

129.0244 

1,99919) 

240 

,4719942 

* , 2505379-02 

.194*372 

344.0*3 

1,714996 

100 

,4094911 

,0444*97 

.2110742 

3*0.9991 

1,432946 

120 

.92749*9 

9.7931744-02 

.2404*11 

377,3393 

1.949*11 

140 

.472*114 

4. 4911999-02 

.2320244 

394.4*17 

2,0*3321 

240 

.4121077 

4.213*739-02 

.2340137 

411,9232 

2.17)04* 

240 

.19737*9 

1.4910499-02 

.2034994 

439.4142 

2,27913 

400 

.  3090144 

3. 2414149-02 

.  1492144 

447,313* 

2,379113 

420 

.2470919 

.0391011 

, 1471223 

4*9.1991 

2,479412 

440 

.2311299 

2.4449429-02 

,131293 

441.3219 

2.9*7413 

440 

.2004419 

3.4949199-03 

. 1349943 

901.3974 

2, *39341 

440 

.1749714 

.0324749 

.  12437 

330 

2,739442 

soo 

,1931799 

2 . 1907149-02 

.1132719 

939.91*9 

3,4199*3 

T«bl«  6.  An  axampla  of  a  lataral  ranga  function  valuaa  printout. 
Tha  haading  for  tha  croaacorralation  valuaa  ia  p(cc) 
and  tha  haading  for  tha  aquara  law  valuaa  ia  p(al) . 
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t  .31*  o'  Ivorvj* 


Figure  5.  A  plot  of  the  square  law  lateral  range  function  values 
that  are  listed  in  Table  6.  The  horizontal  axis  is 
encounter  horizontal  range  at  CPA  in  meters.  The 
vertical  axis  is  encounter  detection  probability. 
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4nttl.Mil  nagnteu?  aignal  valuta 
for  a  lataral  rang*  or  40  matara 


ralativa  CPA 

a  in  aattra 

•87*. *324 
••ii.oaii 

-745,383* 

•701.8488 

-*87,4*9) 

••11.8827 

•870.00* 

•824 . 1884 
•482.2128 
•438.446} 
•3*4,418* 
-jso,  ?->a* 
•306.4263 
•, 47,0  7*7 

•} IV, 2131 

-I'fTi.amtta. 

-131.91*8 

-»7 , k*)24 

*43.846*} 

0 

47,4446} 

67,41)14 

1 U  *  199 
V  "  9  ,  3  3  6  *, 

1 3 , ;  j  )  i 
is) , o’*  > 

3  06, n:i . 

3  90,  '  " 4 
3*4 ,6186 
438,46-} 

49}  .  31.' i 
936 . 19V  * 
970,006 
413,89;' 

6b  7 ,6  99 l 
7  01 , b  4  S  v 
7  4  5  ,  J  9  J 
789,2391 
$33,0858 
876,4314 


Table  7.  An 


diatam.'*  ■agnatic  aignal 

Ha  In  gam* 

4.14*2431-03 
8 .  •88I1S8-03 
1  •  1812411*02 
1 , 8413348-03 
2.0784871-03 
2.7*18788-02 
2.7*72*21-02 
». 1*88771-02 
7.2311*31-02 
.10201*4 
.1444*87 
.2139143 
,317944 
.410489 
.739*84 
1.130141 
1.707447 
2.448238 
3.133584 
3.343334 
3.439438 
1.033839 
-.14*7833 
-.7049798 
-.1055)38 
-.7078842 
-.96314*2 
-.4)07948 
-.335*745 
-.3445983 
-.1881302 
-.1451051 
-.1133414 
-8 . 844089E-03 
-7 . 1929345-02 
-5 . 7122795-02 
-4 .7536425-02 
-3 , 896839E-02 
-3 . 2441735-02 
-2.72)7335-02 
-2,3047615-02 
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-676.9224 


-438.466.;? 


o 


-436.4662 


-676.9324 


Figure  6 


.  A  plot  of  the  magnetic  signal  values  listed  in  Table 
7.  The  horizontal  axis  is  relative  CPA  distance  in 
meters.  The  vertical  axis  is  magnetic  signal  in 


gamma. 


f 

I 


Apptndix  10.  A  Program  Listing 

10  CIS  :  CHAN 

20  PRINT  "Magnetic  Anomaly  Dataction  (MAD)  Lataral  Nanga  Function  Program11 
30  DIM  X0(70),  R0(70),  H(70),  HS(70,  ISO),  K(70),  PDCC(70),  PDSK70) 

AO  PI  •  3. HI 592654#:  CON  ■  2  *  PI  /  360:  KON  "  1B52  /  3600:  N«  ■  "MAD. BAS" 

50  QO  -  .2316419:  91  •  .31938153#:  92  ■  -.356563782#:  93  ■  1.781477937# 

60  94  ■  -1.821255978#:  95  ■  1.330274429# 

70  II  ■  2.515517:  12  -  .802853:  13  ■  .010328:  14  •  1.432788:  15  ■  .189269:  16  ■  .001308 
80  PNINT 

90  INPUT  "ganarata  data  or  print  a  program  data  file  (g/p)";  AS 

100  IF  AS  ■  "G"  ON  AS  ■  "g"  THEN  GOTO  120 

110  IF  AS  ■  "P"  ON  AS  ■  "p"  THEN  GOTO  3070  ELSE  GOTO  90 

120  PNINT  :  AS  »  "a" 

130  INPUT  "magnatic,  procatsing  t  kinematic  data  antry  by  combined  file  (y/n)H;  AS 

140  IF  AS  •  "Y"  ON  AS  ■  "y"  THEN  GOTO  2120 

150  IF  AS  a  «N»  ON  AS  a  Nnii  THEN  qoto  160  ELSE  GOTO  130 

160  INPUT  "magnatic  data  antry  by  file/keyboard  (f/k)";  AS 

170  IF  AS  ■  "F"  ON  AS  ■  "f"  THEN  GOTO  780 

180  IF  AS  ■  "K"  ON  AS  >  a|(ii  THEN  GOTO  190  ELSE  GOTO  160 

190  INPUT  "latituda  in  dacimal  dagraas  (N  ♦>";  LAT:  LAIN  ■  LAT  *  CON 

200  INPUT  "longitude  in  dacimal  degrees  (E  ♦)";  LNG:  LNGR  -  -LNG  *  CON 

210  INPUT  "variation  in  decimal  degrees  (E  ♦  )»;  DEC:  DECR  «  DEC  *  CON 

220  AS  >  »a" 

230  INPUT  "input  dip  angle  (y/n)";  AS 

240  IF  AS  «  »Y»  OR  AS  .  "y"  THEN  GOTO  260 

250  IF  AS  a  »n"  OR  AS  a  "n»  THEN  GOTO  290  ELSE  GOTO  230 

260  INPUT  "dip  angle  in  decimal  degrees  (north  magnetic  hemisphere  ♦)";  DIP 

270  D1PR  a  DIP  •  CON 

280  GOTO  390 

290  LATPR  a  76  *  CON:  LNGPR  -  100  *  CON 

300  D  «  S I N ( LNGR  -  LNGPR)  *  COS(LATR):  E  a  COSCLNGR  -  LNGPR)  *  COS(LATR):  F  *  SIN(LATR) 

310  X  «  E:  Y  a  f 
320  GOSUB  4170 

330  T  a  T  -  (90  *  CON  -  LATPR):  E  *  R  *  SIN(T):  F  •  R  *  COS(T) 

340  X  «  E:  Y  a  d 

350  GOSUB  4170 

360  X  a  F:  Y  a  r 

370  GOSUB  4170 

380  DIPR  a  ATN(2  *  TAN(T)):  DIP  =  D1PR  /  CON 
390  AS  ■  "a" 

400  INPUT  "input  encounter  magnetic  field  intensity  (y/n)";  AS 

410  IF  AS  -  »Y»  OR  AS  «  »y"  THEN  GOTO  430 

420  IF  AS  a  «N"  OR  AS  a  ><n»  THEN  GOTO  450  ELSE  GOTO  400 

430  INPUT  "encounter  magnetic  field  intensity  in  gamma";  HE 

440  GOTO  460 

450  HE  «  70000!  /  SOR(3  *  COS(DIPR)  *  COS(DIPR)  +  1) 

460  PLM  «  0:  PTM  a  0:  PVM  «  0 
470  AS  «  "a" 

480  INPUT  "input  target  permanent  dipole  moments  <y/n)";  AS 

490  IF  AS  a  »y»  OR  AS  *  "y"  THEN  GOTO  510 

500  IF  AS  «  "N"  OR  AS  *  "n"  THEN  GOTO  540  ELSE  GOTO  480 

510  INPUT  "permanent  longitudinal  moment  in  oersted-cm3  (stern-to-bow  ♦)";  PLM  520  INPUT  "permanent 

transverse  moment  in  oersted-cm3  (port-to-starboard  ♦)";  PTM 

530  INPUT  "permanent  vertical  moment  in  oersted-cm3  (downward  ♦)";  PVM 

540  INPUT  "target  displacement  in  tons";  WT 

550  AS  a  "a" 

560  INPUT  "input  target  permeability  coefficients  or  factors  (c/f)";  AS 

570  IF  AS  a  "C»  OR  AS  «  "c"  THEN  GOTO  590 

580  IF  AS  ■  »F"  OR  AS  a  «f«  THEN  GOTO  630  ELSE  GOTO  560 

590  INPUT  "longitudinal  permeability  coefficient  in  cgs  units";  KL 

600  INPUT  "transverse  permeability  coefficient  in  cgs  units";  KT 

610  INPUT  "vertical  permeability  coefficient  in  cgs  units";  KV 

620  FL  «  KL  /  WT:  FT  «  KT  /  WT:  FV  ■  KV  /  WT:  GOTO  710 

630  INPUT  "longitudinal  displacement  factor  in  cgs  units";  FL 
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640  INPUT  "transveraa  displacement  factor  In  cgs  units";  FT 
6S0  INPUT  "varttcat  displacement  factor  In  cga  unlti";  FV 
660  KL  ■  FL  *  WT:  KT  ■  FT  *  WT:  KV  ■  FV  *  WT 
670  AS  ■  "a" 

680  INPUT  "ganarata  a  magnatlc  data  f 1 l a  (y/n)";  AS 

690  IF  AS  ■  "Y"  OR  AS  ■  «y«  THEN  GOTO  710 

700  IF  AS  ■  »N»  Oft  AS  ■  "n"  THEN  GOTO  850  ELSE  GOTO  680 

710  INPUT  "magnetic  data  fila  name";  MS 

720  ON  ERROR  GOTO  730s  GOTO  740 

730  RESUME  710 

740  OPEN  "O",  #1,  MS 

750  WRITE  #1,  LAT,  LNG,  DEC,  DIP,  HE,  PLM,  PTM,  PVH,  WT,  KL,  KT,  KV,  FL,  FT,  FV 
760  CLOSE 
770  GOTO  850 

780  INPUT  "magnatlc  data  fila  nama";  MS 
790  ON  ERROR  GOTO  800:  GOTO  810 
800  RESUME  780 
810  OPEN  "I",  #1,  MS 

820  INPUT  #1,  LAT,  LNG,  DEC,  DIP,  HE,  PLM,  PTM,  PVM,  WT,  KL,  KT,  KV,  FL,  FT,  FV 
830  CLOSE 

840  DECR  ■  DEC  •  CON:  DIPR  ■  DIP  *  CON 
850  AS  ■  "a" 

860  INPUT  "processing  data  entry  by  file/keyboard  (f/k)";  AS 

870  IF  AS  ■  "F"  OR  AS  »  »f"  THEN  GOTO  1450 

880  IF  AS  >  "K"  OR  AS  *  "k"  THEN  GOTO  890  ELSE  GOTO  860 

890  AS  «  "a" 

900  INPUT  "input  sampling  period  (y/n)";  AS 

910  IF  AS  *  «Y"  OR  AS  *  "y"  THEN  GOTO  930 

920  IF  AS  *  »N"  OR  AS  ■  "n"  THEN  GOTO  960  ELSE  GOTO  900 

930  INPUT  "sampling  period  in  seconds";  DT 

940  IF  DT  <■  0  THEN  PRINT  :  PRINT  "must  be  greater  than  zero":  PRINT  :  GOTO  930 
950  GOTO  1090 
960  AS  ■  "a" 

970  INPUT  "input  maximum  magnetic  signal  frequency  (y/n)";  AS 
980  IF  AS  «  "Y»  OR  AS  «  "y"  THEN  GOTO  1000 

990  IF  AS  *  "N"  OR  AS  ■=  "n"  THEN  GOTO  1030  ELSE  GOTO  970 

1000  INPUT  "maximum  magnetic  signal  frequency  in  Hertz";  MAXF 

1010  IF  MAXF  «  0  THEN  PRINT  :  PRINT  "must  be  greater  than  zero":  PRINT  :  GOTO  1000 
1020  GOTO  1080 

1030  INPUT  "minimus  target  slant  range  at  CPA  in  meters";  MINRO 

1040  IF  MINRO  <*  0  THEN  PRINT  :  PRINT  "must  be  greater  than  zero":  PRINT  :  GOTO  1030 
1050  INPUT  "maximun  magnetometer  relative  speed  in  knots";  MAXVMK 
1060  MAXVM  -  MAXVMK  *  KON:  MAXF  ■  2  *  MAXVM  /  MINRO 

1070  IF  MAXVM  <*  0  THEN  PRINT  :  PRINT  "must  be  greater  than  zero":  PRINT  :  GOTO  1050 

1080  DT  *  1  /  (2  *  MAXF):  REM  low  pass  filter,  Nyquist  sampling  rate 

1090  AS  «=  "a" 

1100  INPUT  "input  integration  time  (y/n)";  AS 

1110  IF  AS  ■  »Y"  OR  AS  *  "y"  THEN  GOTO  1130 

1120  IF  AS  *  "N»  OR  AS  *  "n"  THEN  GOTO  1210  ELSE  GOTO  1100 

1130  INPUT  "integration  tire  in  seconds";  IT 

1140  IF  IT  >=  DT  THEN  GOTO  1170 

1150  PRINT  :  PRINT  "IT  *  »  ♦  STRS(IT)  +  »  seconds  minimur  =  »  ♦  STRS(DT)  +  "  seconds":  PRINT 
1160  GOTO  1130 

1170  NS  ■  2  *  INTOT  /  DT  /  2)  ♦  1:  REM  adj  nunber  of  samples  per  integration  time 

1180  IF  NS  <«  151  THEN  GOTO  1340 

1190  PRINT  :  PRINT  "IT  «  »  ♦  STRS(IT)  ♦  "  seconds  maximum  »  »  ♦  STR$(150  *  DT)  ♦  "  seconds":  PRINT 

1200  GOTO  1130 

1210  INPUT  "maximum  target  slant  range  at  CPA  in  meters";  MAXRO 

1220  IF  MAXRO  «  0  THEN  PRINT  :  PRINT  "must  be  greater  than  zero";  PRINT  :  GOTO  1210 
1230  INPUT  "miniaua  magnetometer  relative  speed  in  knots";  MINVMK 
1240  MINVM  -  MINVMK  *  KON 

1250  IF  MINVM  <»  0  THEN  PRINT  :  PRINT  "must  be  greater  than  zero":  FRINT  :  GOTO  1230 
1260  IT  *  2  *  MAXRO  /  MINVM 
1270  IF  IT  »*  DT  THEN  GOTO  1300 

1280  PRINT  :  PRINT  "IT  =  »  ♦  STRS(IT)  ♦  "  seconds  mininur  «  "  ♦  STRS(DT)  ♦  »  seconds":  PRINT 
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1290  GOTO  1100 

1300  NS  ■  2  *  1NT<1T  /  OT  /  2)  +  Is  REM  adjusted  number  of  samples  per  integration  tint 

1310  IF  NS  <>  151  THEN  GOTO  1340 

1320  PRINT  J  PRINT  "IT  •  "  ♦  STRS(IT)  ♦  «  seconds  maxim*  ■  "  ♦  STRSO50  *  DT)  ♦  "  seconds":  PRINT 
1330  GOTO  1100 
1340  AS  ■  "a" 

1350  INPUT  "generate  a  procsssing  data  file  (y/n)";  AS 

1360  IF  AS  >  "Y"  OR  AS  ■  "y"  THIN  GOTO  1380 

1370  IF  AS  ■  «N"  OR  AS  ■  "n"  THIN  GOTO  1510  ELSE  GOTO  1350 

1380  INPUT  "processing  data  file  name";  PS 

1390  ON  ERROR  GOTO  1400:  GOTO  1410 

1400  RESUME  1380 

1410  OPEN  "0",  *1,  PS 

1420  WRITE  *1,  DT,  IT,  NS 

1430  CLOSE 

1440  GOTO  1510 

1450  INPUT  "processing  data  file  name";  PS 

1460  ON  ERROR  GOTO  1470:  GOTO  1480 

1470  RESUME  1450 

1480  OPEN  "I",  «1,  PS 

1490  INPUT  #1,  DT,  IT,  NS 

1500  CLOSE 

1510  AS  ■  "a" 

1520  INPUT  "kinematic  data  entry  by  file/keyboard  (f/k)“;  AS 

1530  IF  AS  ■  "F"  OR  AS  *  "f"  THEN  GOTO  1790 

1540  IF  AS  *  "K"  OR  AS  -  "k"  THEN  GOTO  1550  ELSE  GOTO  1510 


1550  INPUT  "magnetometer  course  in  decimal  degrees  (0  if  at  rest)";  CM 

1560  INPUT  "magnetometer  speed  in  knots";  VMK 

1570  INPUT  "magnetometer  altitude  in  meters  (below  0  is  -)";  AM 

1580  INPUT  "target  course  in  decimal  degrees";  CT 

1590  INPUT  "target  speed  in  knots";  VTK 

1600  INPUT  "target  depth  in  meters  (above  0  is  ->";  AT 

1610  CNR  •  CM  *  CON:  CTR  «  CT  *  CON 

1620  WXK  ■  VMK  *  SIN(CMR)  *  VTK  •  SIN(CTR) 

1630  WYK  ■  VMK  *  COS(CMR)  •  VTK  *  COS(CTR) 

1640  Z  ■  AM  ♦  AT:  REM  vertical  separation  (-  for  magnetometer  below  target) 

1650  X  «  WXK:  Y  «  WYK:  GOSUB  4170 

1660  CR  ■  T:  C  ■  CR  /  CON:  REM  magnetometer  relative  course 
1670  WK  *  R:  REN  magnetometer  relative  speed 
1680  AS  ■  "a" 

1690  INPUT  "generate  a  kinematic  data  file  (y/n)";  AS 

1700  IF  AS  ■  "Y"  OR  AS  »  »y"  THEN  GOTO  1720 

1710  IF  AS  ■  "N"  OR  AS  ■  »n"  THEN  GOTO  1860  ELSE  GOTO  1690 

1720  INPUT  "kinematic  data  file  name";  KS 

1730  ON  ERROR  GOTO  1740:  GOTO  1750 

1740  RESUME  1720 

1750  OPEN  "0",  HU,  KS 

1760  WRITE  81,  CM.  VMK,  AM,  CT,  VTK,  AT,  Z,  C,  WK 
1770  CLOSE 
1780  GOTO  1860 

1790  INPUT  "kinematic  data  file  name";  KS 
1800  ON  ERROR  GOTO  1810:  GOTO  1820 
1810  RESUME  1790 
1820  OPEN  "I",  #1,  KS 

1830  INPUT  81,  CM,  VMK,  AM,  CT,  VTK,  AT,  Z,  C,  WK 
1840  CLOSE 

1850  CMR  ■  CM  *  CON:  CTR  ■  CT  *  CON:  CR  «  C  *  CON 

1860  MCMR  ■  (CMR  -  DECR):  REM  magnetometer  awgnetic  course  in  radians 

1870  MCTR  ■  (CTR  -  DECR):  REM  target  magnetic  course  in  radiana 

1880  ILM  «  KL  *  Ml  •  COS(DIPR)  *  COS(MCTR):  ITM  ■  -KT  *  HE  •  COS(DIPR)  *  SIN(MCTR) 

1890  IVM  «  KV  *  HE  *  SIN(DIPR) 

1900  DMX  ■  (PLM  ♦  ILM)  *  SIN(MCTR)  ♦  (PTM  +  ITM)  •  COS(NCTR) 

1910  DMY  -  (PLM  ♦  ILM)  *  COS(MCTR)  -  (PTM  ♦  ITM)  •  SIN(MCTR) 

1920  DMV  «  PVM  ♦  IVM 
1930  X  *  DMX:  Y  ■  DMY 
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1940  GOSUB  4170 

1950  OMLR  •  T:  REM  dipole  azimuth  rtlatlv*  to  magnetic  north 
1960  X  -  OMV:  Y  «  R 
1970  GOSUB  4170 

1960  DM  *  R:  OMR  ■  T:  REM  dipole  depreaaion  angle 
1990  OML  *  OMLR  /  CON:  OM  ■  OMR  /  CON 
2000  AS  ■  "a" 

2010  INPUT  "generate  e  combined  magnetic,  proceaalng  t  kinematic  data  file  (y/n)H;  AS 

2020  IF  AS  ■  "Y"  OR  AS  ■  "y"  THEN  GOTO  2040 

2030  IF  AS  ■  "N"  OR  AS  ■  "n"  THEN  GOTO  2210  ELSE  GOTO  2010 

2040  INPUT  "combined  magnetic,  proceaalng  t  kinematic  data  file  name";  ES 

2050  ON  ERROR  GOTO  2060:  GOTO  2070 

2060  RESUME  2040 

2070  OPEN  "0",  #1,  ES 

2080  WRITE  #1,  LAT,  LNG,  DEC,  DIP,  HE,  PLM,  PTM,  PVM,  UT,  KL,  KT,  KV,  FL,  FT,  FV,  DT,  IT 
2090  WRITE  #1,  NS,  CM,  VMK,  AM,  CT,  VTK,  AT,  2,  C,  WK,  ILM,  ITM,  IVM,  DM,  OML,  OM,  MS,  PS,  ICS 
2100  CLOSE 
2110  GOTO  2210 

2120  INPUT  "combined  magnetic,  proceaalng  t  kinematic  data  file  name";  ES 
2130  ON  ERROR  GOTO  2140:  GOTO  2150 
2140  RESUME  2120 
2150  OPEN  "I",  #1,  ES 

2160  INPUT  #1,  LAT,  LNG,  DEC,  DIP,  HE,  PLM,  PTM,  PVM,  WT,  KL,  KT,  KV,  FL,  FT,  FV,  DT,  IT 
2170  INPUT  #1,  NS,  CM,  VMK,  AM,  CT,  VTK,  AT,  2,  C,  WK,  ILM,  ITM,  IVM,  DM,  OML,  OM,  MS,  PS,  KS 
2180  CLOSE 

2190  DECR  «  DEC  *  CON:  OIPR  *  DIP  *  CON 

2200  CR  «  C  *  CON:  OMLR  ■  OML  *  CON:  OMR  *  OM  *  CON 

2210  INPUT  "required  falae  alarm  rate  in  false  alarms  per  hour";  FAR 

2220  PF  ■  FAR  *  IT  /  3600:  REM  falae  alarm  probability 

2230  Y  ■  PF:  IF  PF  >  .5  THEN  Y  ■  1  -  Y:  REM  inverse  normal  approximation 

2240  Y  ■  SORCLOCCI  /  Y  /  Y>) 

2250  Y  ■  Y  -  (II  ♦  Y  *  <12  ♦  13  *  Y>)  /  <1  ♦  Y  *  <14  ♦  Y  *  (15  ♦  16  *  Y)>) 

2260  IF  PF  <  .5  THEN  Y  «  -Y 

2270  ZP  ■  -Y 

2280  CHI  ■  NS  *  (1  •  2  /  9  /  NS  ♦  ZP  *  SQR(2  /  9  /  NS))  *  3:  REM  inverse  chi-square  approximation 

2290  INPUT  "magnetic  noise  in  gamna";  SIG 

2300  INPUT  "maximum  lateral  range  in  meters";  LRM 

2310  INPUT  "lateral  range  step  in  meters";  ST 

2320  IF  ST  <»  LRM  THEN  GOTO  2340 

2330  PRINT  :  PRINT  “maximum  step  is  "  ♦  STRS(LRM)  ♦  »  meters":  PRINT  :  GOTO  2310 
2340  NC  ■  2  *  INTCLRM  /  ST)  ♦  1:  REM  nunber  of  lateral  range  function  values 
2350  IF  NC  «  71  THEN  GOTO  2370 

2360  PRINT  :  PRINT  "minimum  step  is  "  ♦  STRSCLRM  /  35)  ♦  "  meters":  PRINT  :  GOTO  2310 

2370  AIT  ■  DT  *  NS:  REM  adjusted  integration  time 

2380  W  ■  WK  *  KON:  REM  magnetometer  relative  speed  in  meters/second 

2390  DS  ■  W  *  DT:  REM  distance  between  samples  on  the  relative  track 

2400  XO  -  -(NC  -  1)  /  2  •  ST 

2410  FOR  I  «  0  TO  NC  -  1 

2420  X0( I >  *  XO 

2430  X  «  XO:  Y  ■  2 

2440  GOSUB  4170 

2450  RO  *  R :  R0( I )  ■  R:  REM  target  slant  range  at  CPA  in  meters 
2460  DELR  *  T:  REM  target  depression  angle  complement  at  CPA  in  radians 

2470  IF  RO  •  0  THEN  GOTO  2800:  REM  zero  lateral  range  and  vertical  separation  2480  DMF  ■  DM  /  10  /  RO  * 
REM  dipole  moment  factor 

2490  CMR  ■  CR  -  OFCR:  REM  target  relative  magnetic  course 
2500  L  ■  COS(OMR)  *  COS (CMR  -  OMLR) 

2510  M  >  COS(DELR)  *  COS(OMR)  *  SINCCMR  -  OMLR)  -  SIN(DELR)  *  SIN(OMR) 

2520  N  ■  -SIN(OELR)  *  COS(OMR)  •  SINCCMR  -  OMLR)  -  COS(DELR)  *  SIN(OMR) 

2530  LI  «  COS(DIPR)  *  COS(CMR) 

2540  Ml  «  COS(OELR)  *  COS(OIPR)  *  SINCCMR)  •  SIN(DELR)  *  SIN(DIPR) 

2550  N1  ■  -SIN(DELR)  *  COS(DIPR)  *  SINCCMR)  •  COSCDELR)  *  S1NCDIPR) 

2560  A2  *  2  *  L  *  Li  -  M  *  Ml  -  N  *  N1:  REM  Anderson  Function  Coefficient 
2570  A1  «  3  *  (N  *  LI  ♦  L  *  Nl):  REM  Anderson  Function  Coefficient 
2580  AO  •  2  *  N  *  Nl  -  L  *  LI  -  M  *  Ml:  REM  Anderson  Function  Coefficient 
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2590  SUN  ■  Os  HNAX  ■  0:  NMtN  ■  0 
2600  FOR  J  ■  0  TO  NS  -  1 

2610  S  ■  <J  *  (MS  •  1)  /  2)  *  OS:  BA  «  S  /  RO:  REN  Anderson  Function  Argument  2620  AF  -  1  /  (1  ♦  BA  *  BA)  A 


2.5:  REM  Anderaon  Function  Factor 

2630  NSF  ■  <A2  *  BA  *  BA  +  A1  •  BA  ♦  AO)  •  AF:  REM  magnetic  signal  factor 
2640  HSU,  J)  ■  DNF  *  NSF:  REM  magnetic  signal  value 
2650  IF  HS( I ,  J)  >  HNAX  THEN  HNAX  >  HS(I,  J) 

2660  IF  HSU,  J)  <  HMIM  THEN  HM1N  «  HSU,  J) 

2670  SUN  ■  SUM  ♦  HSF  *  HSF 

2680  NEXT  J 

2690  HU)  ■  HMAX  -  HMIN 
2700  KU)  ■  SOR(SUN) 

2710  W  •  -2P  ♦  DNF  *  SOR(SUM)  /  S1G 

2720  LAN  ■  ONF  *  DNF  *  SUN  /  (S1G  *  SIG):  AA  ■  NS  +  LAM:  BB  «  1  ♦  LAN  /  (NS  «•  LAM) 

2730  ZN  »  -SQR(2  *  CHI  /  BB)  *  SOR(2  *  AA  /  BB  -  1):  XI  »  W 

2740  GOSUB  4230 

2750  IF  Y1  >  1  THEN  Y1  ■  1 

2760  POCCU)  -  Y1 :  XI  ■  ZN 

2770  GOSUB  4230 

2780  IF  Y1  >  1  THEN  Y1  *  1 

2790  PDSLU)  ■  Y1 

2800  XO  «  XO  «■  ST 

2810  NEXT  I 

2820  C  «  C  /  360:  C  *  (C  -  INT(C))  *  360 
2830  IF  C  <  0  THEN  C  =  360  ♦  C 

2840  OML  -  (OML  +  DEC)  /  360:  OML  •  (OML  •  INT(OML))  *  360 
2850  IF  OML  <  0  THEN  OML  *  360  ♦  OML 
2860  PRINT  :  AS  ■  "a" 


2870 

2880 

2890 

2900 

2910 

2920 

2930 

2940 

2950 

2960 

2970 

2980 

2990 

3000 

3010 

3020 

3030 

3040 

3050 

3060 


INPUT  "generate  a  program  data  file  (y/n)H;  AS 

IF  AS  «  «Y"  OR  AS  «  "y"  THEN  GOTO  2900 

IF  AS  ■  "N"  OR  AS  «  "n"  THEN  GOTO  3230  ELSE  GOTO  2870 

INPUT  "program  data  file  name";  DS 

ON  ERROR  GOTO  2920:  GOTO  2930 

RESUME  2900 

OPEN  "0",  #1,  DS 

WRITE  #1,  LAT,  LNG,  DEC,  DIP,  HE,  PLM,  PTM,  PVN,  WT,  XL,  XT, 
NS,  CM,  VMX,  AM,  CT,  VTX,  AT,  ILN,  ITN,  IVM,  Z,  C, 
LRM,  DS,  NC,  MS,  PS,  XS,  ES 
TO  NC  -  1 

X0(1),  POCCU),  PDSLU),  XU),  HU),  R0( I ) 


WRITE  #1, 
WRITE  #1, 
FOR  I  «  0 
WRITE  *1, 
NEXT  I 
FOR  I  »  0 
FOR  J  «  0 
WRITE  #1, 
NEXT  J 
NEXT  I 
CLOSE 
GOTO  3230 


TO  NC  -  1 
TO  NS  •  1 
HSU,  J) 


3070  INPUT  "program  data  file  name";  DS 
3080  ON  ERROR  GOTO  3090:  GOTO  3100 
3090  RESUME  3070 
3100  OPEN  "I",  #1,  DS 

3110  INPUT  HU,  LAT,  LNG,  DEC,  DIP,  HE,  PLM,  PTM,  PVN,  UT,  XL,  XT, 
3120  INPUT  #1,  NS,  CM,  VMX,  AM,  CT,  VTX,  AT,  ILM,  ITM,  IVM,  Z,  C, 

3130  INPUT  #1,  LRM,  DS,  NC,  MS,  PS,  XS,  ES 

3140  FOR  I  •  0  TO  NC  -  1 

3150  INPUT  #1,  X0(1),  POCCU),  PDSLU),  XU),  HU),  R0{I) 

3160  NEXT  I 

3170  FOR  I  •  0  TO  NC  -  1 

3180  FOR  J  •  0  TO  NS  •  1 

3190  INPUT  #1,  HSU,  J) 

3200  NEXT  J 
3210  NEXT  I 
3220  CLOSE 

3230  PRINT  :  AS  >  "a" 


XV,  FL,  FT,  FV,  DT,  IT,  AIT 
WX,  DM,  OML,  OM,  FAR,  PF,  SIG,  ST 


XV,  FL,  FT,  FV,  DT,  IT,  AIT 
WK,  DM,  OML,  OM,  FAR,  PF,  SIG,  ST 
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3240  INPUT  "print  encounter  ptriMttr  values  (y/n)";  AS 
3250  IP  AS  ■  "Y"  OR  AS  •  "y"  THEN  GOTO  3270 

3260  IP  AS  «  "N»  OR  AS  ■  "n«  THEN  GOTO  3780  ELSE  GOTO  3240 

3270  LPRINT 

3280  LPRINT  "program  file  name  "  ♦  NS 

3290  LPRINT  "program  data  Pi  la  name  "  *  OS 

3300  LPRINT  "Magnetic  data  file  name  "  ♦  NS 

3310  LPRINT  "processing  data  file  name  "  ♦  PS 

3320  LPRINT  "kinemtic  data  file  name  "  ♦  KS 

3330  LPRINT  "combined  magnetic,  processing  t  kinematic  data  file  name  "  ♦  ES 

3340  LPRINT  "encounter  latitude  (decimal  degrees)  ";  SPC(2); 

3350  LPRINT  "encounter  longitude  (decimal  degrees)  ";  SPC(2>; 

3360  LPRINT  "encounter  variation  (decimal  degrees)  ";  SPC(2); 

3370  LPRINT  "encounter  dip  angle  (decimal  degrees)  ";  SPC(2); 

3380  LPRINT  "encounter  magnetic  field  intensity  (gamma)  ";  SPC(2); 

3390  LPRINT  "permanent  longitudinal  moment  (oersted-cm3)  ";  SPC(2); 

3400  LPRINT  "permanent  transverse  moment  (oersted-cm3)  ";  SPC(2); 

3410  LPRINT  "permanent  vertical  moment  (oersted-cm3)  ";  SPC(2); 

3420  LPRINT  "target  displacement  (tons)  ";  SPC(2); 

3430  LPRINT  "target  longitudinal  permeability  coefficient  ";  SPC(2); 

3440  LPRINT  "target  transverse  permeability  coefficient  ";  SPC(2); 

3450  LPRINT  "target  vertical  permeability  coefficient  ";  SPC(2); 

3460  LPRINT  "target  longitudinal  permeability  factor  ";  SPC(2); 

3470  LPRINT  "target  transverse  permeability  factor  SPC(2); 

3480  LPRINT  "target  vertical  permeability  factor  ";  SPC(2); 

3490  LPRINT  "sampling  period  (seconds)  ";  SPC(2); 

3500  LPRINT  "integration  time  (seconds)  SPC(2); 

3510  LPRINT  "adjusted  integration  time  (seconds)  ";  SPC (I); 

3520  LPRINT  "number  of  sanples  per  encounter  SPC(2); 

3530  LPRINT  "magnetometer  course  (decimal  degrees)  ";  SPC(2); 

3540  LPRINT  "magnetometer  speed  (knots)  ";  SPC(2); 

3550  LPRINT  "magnetometer  altitude  (meters)  ";  SPC(2); 

3560  LPRINT  "target  course  (decimal  degrees)  ";  SPC(2); 

3570  LPRINT  "target  speed  (knots)  ";  SPC(2); 

3580  LPRINT  “target  depth  (meters)  »;  SPC(2); 

3590  LPRINT  "magnetometer  relative  course  (decimal  degrees)  SPC(2); 

3600  LPRINT  "magnetometer  relative  speed  (knots)  ";  SPC(2); 

3610  LPRINT  "magnetometer- target  vertical  separation  (meters)  ";  SPC(?>; 

3620  LPRINT  "target  induced  longitudinal  dipole  moment  (oersted-cm3)";  SPC(2j; 

3630  LPRINT  "target  induced  transverse  dipole  moment  (oersted-cm3)  SPC(2); 

3640  LPRINT  "target  induced  vertical  dipole  moment  (oersted-cm3)  SPC(2); 

3650  LPRINT  "magnetic  dipole  moment  (oersted-cm3)  ";  SPC(2); 

3660  LPRINT  "dipole  moment  azimuth  (decimal  degrees)  SPC(2); 

3670  LPRINT  "dipole  moment  depression  angle  (decimal  degrees)  ";  $PC(2); 

3680  LPRINT  "distance  between  samples  on  the  relative  track  (meters)";  SPC(2); 

3690  LPRINT  "false  alarm  rate  (false  alarms  per  hour)  ";  SPC(2); 

3700  LPRINT  "false  alarm  probability  SPC(2); 

3710  LPRINT  "magnetic  noise  (gamna)  SPC(2); 

3720  LPRINT  "maximum  lateral  range  (meters)  SPC(2); 

3730  LPRINT  “lateral  range  step  (meters)  ";  SPC(2); 

3740  LPRINT  "rxjnber  of  lateral  range  function  values  ";  SPC(2); 

3750  FOR  I  «  0  TO  15 
3760  LPRINT 
3770  NEXT  I 
3780  PRINT  :  AS  «  "a" 

3790  INPUT  "print  lateral  range  function  values  (y/n)M;  AS 
3800  IF  AS  «  "Y"  OR  AS  «  "y"  THEN  GOTO  3820 
3810  IF  AS  •  "N"  OR  AS  «  "n"  THEN  GOTO  3900  ELSE  GOTO  3790 
3820  LPRINT  DS;  "  lateral  range  function  values" 

3830  LPRINT  :  LPRINT 

3840  LPRINT  »L  p(cc)  p(sl)  H  R 

3850  LPRINT  "meters  gamma  meters" 

3860  LPRINT 

3870  FOR  I  «  0  TO  NC  -  1 

3880  LPRINT  X0(I);  TAB(10);  P0CC(1);  TA8(24);  POSL(I);  TAB(38>;  H(I);  TAB(52); 


"  ♦  NS 

"  ♦  DS 

"  ♦  NS 

"  ♦  PS 

"  ♦  KS 

*  ♦  ES 

SPC(2); 

LAT 

S?C<2); 

LNG 

SPC(2); 

DEC 

SPC(2); 

DIP 

SPC(2); 

HE 

SPC(2); 

PLN 

SPC(2); 

PTH 

SPC(2); 

PVM 

SPC(2); 

WT 

SPC(2); 

KL 

SPC(2); 

KT 

SPC(2); 

KV 

SPC(2); 

FL 

SPC(2); 

FT 

SPC(2); 

FV 

SPC(2); 

DT 

IPC(2); 

IT 

BPC(I); 

AIT 

SPC(2); 

NS 

SPC(2); 

CM 

SPC(2); 

VMK 

SPC(2); 

AM 

SPC(2); 

CT 

SPC(2); 

VTK 

SPC(2); 

AT 

SPC(2); 

C 

SPC(2); 

UK 

SPC(?>; 

Z 

SPC(2); 

ILM 

SPC(2); 

ITM 

SPC(2); 

IVW 

SPC(2); 

DM 

SPC(2); 

OML 

SPC(2); 

OH 

SPC(2); 

DS 

SPC(2); 

FAR 

SPC(2); 

PF 

SPC(2); 

SIG 

SPC(2); 

LRM 

SPC(2); 

ST 

SPC(2>; 

NC 

R0( I );  TAB(70);  K(I) 
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3890  NEXT  I 

3900  PRINT  :  AS  ■  "a" 

3910  INPUT  "print  magnetic  aignal  values  (y/n)";  AS 

3920  IF  AS  ■  "Y«  OR  AS  -  "y"  THEN  GOTO  3940 

3930  IF  AS  ■  "N"  OR  AS  «  "n"  THEN  GOTO  4120  ELSE  GOTO  3910 

3940  PRINT  :  PRINT  "Identify  signal  by  ancounter  lateral  range  index" 

3950  PRINT  "lateral  range  equals  the  Index  times  "  ♦  STRS(ST)  ♦  "  maters" 

3960  PRINT  "Index  values:  •"  ♦  STRStCNC  -  1)  /  2)  ♦  "  to  "  ♦  STRSUNC  -  1)  /  2) 

3970  INPUT  "lateral  range  Index";  K:  I  «  K  ♦  <NC  -  1)  /  2 

3980  LPRINT  :  LPRINT  :  L PRINT 

3990  LPRINT  DS;  "  magnetic  signal  values" 

4000  LPRINT  "for  a  lateral  range  of  ";  X  *  ST;  "  meters" 

4010  LPRINT 

4020  LPRINT  "relative  CPA  distance";  TAB(35);  "magnetic  signal" 

4030  LPRINT  "s  In  meters";  TAB<35>;  "Hs  in  ganrna" 

4040  LPRINT 

4050  FOR  J  *  0  TO  NS  •  1 

4060  LPRINT  <4  -  (NS  -  1>  /  2)  *  OS;  TAB<35);  HS(I,  J) 

4070  NEXT  J 

4080  PRINT  :  AS  ■  "a" 

4090  INPUT  "print  magnetic  aignal  values  for  a  different  lateral  range  (y/n)";  AS 

4100  IF  AS  ■  "V"  OR  AS  «  "y"  THEN  GOTO  3940 

4110  IF  AS  ■  "N"  OR  AS  ■  "n"  THEN  GOTO  4120  ELSE  GOTO  4090 

4120  PRINT  :  AS  *  "a" 

4130  INPUT  "continue  to  use  the  program  (y/n)";  AS 
4140  IF  AS  «  "Y"  OR  AS  =  "y"  THEN  GOTO  10 

4150  IF  AS  *  "N"  OR  AS  «  "n"  THEN  GOTO  4160  ELSE  GOTO  4130 

4160  ENO 

4170  R  ■  SQR(X  *  X  ♦  Y  •  Y):  REM  rectangular  to  polar  conversion 
4180  IF  R  «  0  THEN  T  >  0:  RETURN 

4190  IF  A8S(X  /  R)  ■  1  THEN  Q  »  SGN(X)  *  (PI  /  2)  ELSE  0  ■  ATN(X  /  R  /  SQR(1  -  X  *  X  /  R  /  R>) 

4200  IF  AB$(Y  /  R)  «  1  THEN  T  ■  (PI  /  2)  *  (1  -  SGN(Y))  ELSE  T  «  (PI  /  2)  -  ATN(Y  /  R  /  SQR(1  -  Y  *  Y  /  R  / 

R» 

4210  IF  0  <  0  THEN  T  *  2  *  PI  -  T 
4220  RETURN 

4230  Y1  «  XI:  IF  XI  «  0  THEN  Y1  «  -VI :  REH  normal  approximation 
4240  G  «  1  /  (1  ♦  OO  *  Y1) 

4250  Yl  ■  EXPt-YI  *  Y1  /  2)  /  SORC2  *  PI )  *  G  *  (Ql  ♦  G  *  (02  ♦  G  *  (03  ♦  G  *  (04  +  G  *  Q5)))) 

4260  IF  XI  >«  0  THEN  Yl  *  1  -  Yl 
4270  RETURN 
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